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We show that in hybrid inflation it is possible to generate large second-order perturbations in 
the cosmic microwave background due to the instability of the tachyonic field during preheating. 
We carefully calculate this effect from the tachyon contribution to the gauge-invariant curvature 
perturbation, clarifying some confusion in the literature concerning nonlocal terms in the tachyon 
curvature perturbation; we show explicitly that such terms are absent. We quantitatively compute 
the nongaussianity generated by the tachyon field during the preheating phase and translate the 
experimental constraints on the nonlinearity parameter Jnl into constraints on the parameters of 
the model. We also show that nonscale-invariant second-order perturbations from the tachyon field 
with spectral index n — 4 can become larger than the inflaton-generated first-order perturbations, 
leading to stronger constraints than those coming from nongaussianity. The width of the excluded 
region in terms of the logarithm of the dimensionless coupling g, grows linearly with the log of the 
ratio of the Planck mass to the tachyon VEV, log(M p /u); hence very large regions are ruled out if 
the inflationary scale v is small. We apply these results to string-theoretic brane-antibrane inflation, 
and find a stringent upper bound on the string coupling, g s < 10 -4 ' 5 . 

PACS numbers: 11.25. Wx, 98.80.Cq 



I. INTRODUCTION 

Inflation has become the dominant paradigm for the 
generation of the anisotropics of the Cosmic Microwave 
Background (CMB) (see 013 for a review). During infla- 
tion quantum fluctuations are redshifted out of the hori- 
zon where they become "frozen" and later re-enter the 
horizon after Big Bang Nucleosynthesis. Observation of 
the temperature anisotropies of the CMB directly con- 
strain the curvature perturbation somewhat before cos- 
mological scales re-enter the horizon. Current observa- 
tions show that the curvature perturbation is gaussian 
with a scale-invariant spectrum to within experimental 
limits, which is consistent with the predictions of infla- 
tion. Nevertheless, some nongaussianity is expected to be 
generated during inflation 3-[8j (see [!j for a review). It 
is interesting to study the level of nongaussianity usually 
parameterized by a dimensionless nonlinearity parameter 
/nl, in various scenarios for inflation, since it can be an 
important way to discriminate between different models. 
For example, the curvaton mechanism [lj} f° r the gener- 
ation of cosmological perturbations may be ruled out by 
future non-observation of nongaussianity (see, how- 
ever, [13). 

One way to generate significant levels of nongaussian- 
ity at the end of inflation is by preheating |l3 | -[l ^ |. 
Preheating refers to the nonperturbative production of 
^articles which arises either due to coherent oscillations 
20] or tachyonic (spinodal) instability 22]- 24] (sec also 
251 l2t|). It has also been suggested that preheating 
might affect the first order variation of the curvature per- 
turbation, perhaps dominating the contribution of the 
inflaton and thus modifying the standard predictions of 
inflation [^-[Ij. 



In this paper we carefully reconsider the generation 
of nongaussian fluctuations from tachyonic preheating in 
hybrid inflation |31| using second order cosmological per- 
turbation theory similarly to [|| l^]. See j3^| for con- 
straints on the parameter space of hybrid inflation com- 
ing from the WMAP data. During hybrid inflation the 
inflaton field if is displaced from its preferred value, and 
inflation is driven by the false vacuum energy of the "wa- 
terfall" or tachyon field a which is trapped in its false 
vacuum by interactions with the inflaton. When the in- 
flaton reaches some critical value the tachyon effective 
mass squared becomes negative and its fluctuations grow 
exponentially due to the spinodal instability. The ex- 
ponential growth continues until the fluctuations of the 
tachyon start to oscillate about the true vaccuum. At this 
stage the tachyon fluctuations become nonperturbative 
and their back-reaction brings inflation to an end. Here 
we study the evolution of the second-order curvature per- 
turbation due to the tachyonic instability and find that it 
can lead to a large level of nongaussianity, if the tachyon 
was lighter than the Hubble scale m a < 3H/2, during a 
sufficient part of the inflationary epoch. (If m a H, 
the fluctuations of a are exponentially suppressed dur- 
ing inflation.) In the hybrid inflation model, it can be 
natural to have m CT < 3H/2 during inflation, since nv\ 
changes sign at the end of inflation. If the inflaton is 
rolling slowly enough, rr? a will naturally already be close 
to zero even before the end of inflation, at the time of 
horizon crossing. 

The scenario is summarized as follows. During infla- 
tion the squared mass of the tachyon field varies linearly 
with the number of e-foldings of inflation, in the vicinity 
of the point where it passes through zero: m 2 — cH 2 N, 
taking N = to be the value where m^. = 0, and c is a 
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function of the parameters of hybrid inflation. Inflation 
starts at some N = iVj < 0, and ends at N = > 0, 
for a total of N e = —Ni + N* e-foldings. We are going 
to limit our inquiry to regions of parameter space where 
the linear behavior of m 2 is valid during the full du- 
ration of inflation, since this technical assumption makes 
the calculations tractable; however we will show that this 
assumption must often be satisfied anyway, as a conse- 
quence of the experimental limit on the inflationary spec- 
tral index. At N = 0, fluctuations of the tachyon field 
start to grow exponentially, and can make a large con- 
tribution to the curvature perturbation at second order 
in cosmological perturbation theory, resulting in a scale- 
noninvariant component both in the CMB temperature 
power spectrum, and its nongaussian correlations. Our 
goal is to derive constraints on the parameter space of 
hybrid inflation from these effects. 

We start in section[n]by defining the perturbations up 
to second order in the metric and inflationary fields. In 
section lTTTl the hybrid inflation model is reviewed, starting 
with the dynamics of the fields at zeroth order, and then 
their first order perturbations, with emphasis on the dy- 
namics of the tachyon field fluctuations toward the end of 
inflation. In section lTVI we solve for the second order cur- 
vature perturbation which is induced by the tachyonic 
fluctuations. This is used in section to compute the 
bispectrum (three-point function) and the nonlinearity 
parameter /ml, as well as contributions to the spectrum 
itself. In section IVII we incorporate experimental con- 
straints on nongaussianity as well as on the spectrum, to 
derive excluded regions in the parameter space of the hy- 
brid inflation model. In section fVlII we consider the pos- 
sibility of generating significant nongaussianity in brane 



inflation. We conclude in section IVIIII comparing our 
results with previous treatments. We give technical de- 
tails concerning mode functions in de Sitter space, the 
second order perturbed Einstein equations, the inflaton 
curvature perturbation, Fourier transforms of convolu- 
tions, the construction of the tachyon curvature pertur- 
bation and the approximation of an adiabatically varying 
tachyon mass in appendices A-F. 



II. METRIC AND MATTER PERTURBATIONS 



In this section we write down the perturbations of the 
metric and matter fields about a spatially flat Robertson- 
Walker background following 6] . Greek indices run over 
the full spacetime fx, v = 0, 1, 2, 3 while latin indices run 
only over the spatial directions i, j — 1,2, 3. We will work 
mostly in conformal time r, related to cosmic time t by 
dt = adr. Differentiation with respect to conformal time 
will be denoted by /' = d T f and with respect to cosmic 
time by / = d t f. 

The metric is expanded up to second order in fluctua- 



tions as 

5oo = ~a{r) 

g i = a(r) 2 

9ij = a(r) 2 
1 
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where D t j = d L dj — ^Sijd k dk is a trace-free operator. The 
fluctuations are decomposed such that the vector pertur- 

(21 ' (2) 

bations are transverse d 1 ^ = d % x\ — while the 
tensor perturbations are transverse, traceless and sym- 

P) - n ,, j (2) _ n .,(2) (2) 



metric: d l x£' = 0, xT' = °> x\]' = X yl 



In the above 

we have neglected the vector and tensor perturbations at 
first order, which are small, since vector perturbations 
decay with time, while tensors are suppressed by the 
slow roll parameter e. The same is not true at second 
order, however, since the second order tensors and vec- 
tors are sourced by the first order scalar perturbations. 
We adopt the generalized longitudinal gauge defined by 



,(2) 



(2) 



and x (1) = x {2) = 0. The metric 



in this gauge becomes 



9oa 
9oi 
9ij 



-a{rf 



1 + 20 
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= a(r) 2 [(1 - 2V> (1) - ^ (2) ) S k 



+ -(dixf+dixT'+x)?) 



(2) , „(2)> 



(4) 
(5) 
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In hybrid inflation the matter content consists of two 
scalar fields which are expanded in perturbation theory 
as 

ip(r,x) = <p (t) + 5 {1) <p{t,x) + ^5 {2) ip(T,x) (7) 
<j(t,x) - a Q (T)+S (1 ^a(T,x) + ^a{T,x). (8) 

where <p is the inflaton and a the tachyon (or "waterfall" 
field). In hybrid inflation the time-dependent vacuum 
expectation value (VEV) of the tachyon field is set to 
zero co(t) = 0, about which we will say more later. 

The perturbations are defined so that (6®<p) = 0, 
hence (<^(r,x)) = ipo(r). At first order in perturbation 
theory this is automatic since 5^ if contains only one an- 
nihilation/creation operator. However, at higher order in 
perturbation theory the homogeneous k = mode of the 
fluctuation must be subtracted by hand in order to en- 
sure that all of the zero mode of the field is described by 
the nonperturbation background. 

The Einstein tensor and stress-energy tensor expanded 
up to second order in perturbation theory can be found 
in [33|. We do not reproduce these results here, but we 
have carefully checked all the results from which are 
relevant for our analysis. 
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III. HYBRID INFLATION 

We consider hybrid inflation in which both the inflaton 
and the tachyon are real fields with the potential 



V { l p,o) = \^-v>f + ^ 



<£_ 

2 



^V. (9) 



This potential will give rise to topological defects at 
the end of inflation — domain walls in the a field — which 
could produce large nongaussianities apart from the ones 
which we consider (Hi))- But these domain walls must 
somehow be rendered unstable to avoid overclosure of 
the universe, for example through the addition of a small 
term like n<7 3 to JJjJ. We will ignore this issue here. It 
would be interesting to consider a complex tachyon field 
[(57j whose defects formed at the end of inflation are cos- 
mic strings, which are much more phenomenologically 
viable than domain walls. 



We will make extensive use of the slow roll parameters, 
defined by 



e - r) 



so that, during inflation 
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T] = 4 



Xv 4 • 



M p m 2 
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(13) 



Notice that if to 2 tp 2 , <C Aw 4 , then e <C r\- This is equiv- 
alent to demanding that the false vacuum energy of the 
tachyon dominates during inflation, which is the assump- 
tion usually made for hybrid inflation: 



A. Background Dynamics 

At the homogeneous level, the usual Friedmann and 
Klein-Gordon equations for the scale factor and the mat- 
ter fields are 



Xv 4 



Xv 4 



3H 2 = _(^ + ag)+« a V, 
dV 

= cpo + 'SHtpQ + — 
dip 

dV 

= a Q +3H& Q + -5- 



(10) 

(11) 
(12) 



where k 2 = Af~ 2 = 8itGn- Here and elsewhere the po- 
tential and its derivatives are understood to be evalu- 
ated on the background values of the fields so that V — 
V((p ,a ) and dV/dtp = dV/d(p\ {tp=tpo ^ =tTo} , for exam- 
ple. For the potential © we have dV/da = d 2 V / dad<p = 
0, provided that do = 0. We will apply this simplification 
to all subsequent results. 

To see why one should set <7o(t) = 0, notice that the 
tachyon effective mass is m 2 = g 2 p\ — Xv 2 > 0. That 
is, the tachyon mass squared starts out being positive 
during inflation. Provided that there was a long enough 
prior period of inflation, any initial departure of a from 
zero would be exponentially damped. At some point, 
m 2 becomes negative, and the tachyonic instability be- 
gins. However, it is still true that co(t) remains zero 
even then, since the universe will consist of equal num- 
bers of domains with a > and a < 0. On average, these 
give zero, which is the definition of the zeroth order field 
(To • The departures of a from zero between domain walls 
which form are taken account in the fluctuations of the 
field. Thus it is consistent to set <tq(t) = for all times 
in our analysis. (Notice that the situation would be dif- 
ferent if we considered a complex scalar field a in which 
case the phase transition would lead to cosmic string for- 
mation and the radial degree of freedom \a\ would not 
average to zero at late times.) 



^ 0l<T0 = 0) = ^ r + -^£^. 
During the slow roll phase the inflaton equation of mo- 



tion 3H<po + m 2 tpQ = has solution 



ip (t) = </? s exp 



m 2 v (t-t s ) 
3H 



<Ps 



^V" (14) 
a* J 



where we used a(t) = a s e H ^~ ts \ with t s an arbitrary 
time. The Hubble scale remains approximately constant, 
3-ff 2 = Xv 4 /(AM 2 ). Since (fo is decreasing and ao = 0, 
the slow roll parameter e actually decreases slowly during 
inflation while r\ remains constant. 



B. Inflationary Dynamics of First Order 
Fluctuations 



Having argued that the background inflaton field ipo is 
described by l|14f> with H approximately constant during 
both inflation and the tachyonic instability phase, we now 
briefly discuss the dynamics of the first order metric and 
inflaton fluctuations for e,ry <C 1. This section is largely 
review since when cto = the first order metric and infla- 
ton perturbations obey exactly the same equations as in 
single field inflation; however, we include some details of 
the calculations since similar equations will arise when we 
study the second order metric fluctuations. 1 We work in 
conformal time which is most convenient for the ensuing 
calculations. 



1 As will be shown, the second order metric fluctuation cj>( 2 ' obeys 
an inhomogeneous equation where the differential operator is 
identical to the one which determines the dynamics of <f>^; thus 
an understanding of the first order solutions simplifies the con- 
struction of the Green function for the second order fluctuations. 
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The 5^G) = K 2 S^Ti Einstein equation for i ^ j 

implies that (j^ 1 ' = ip^ 1 ' which is a well known result. We 
apply this simplification in the following. The S^G® = 
K 2 j(i) j-o Einstein equation is a constraint 



+ H4> 



(i) 



(15) 



which means that the first order metric perturbation (M > 
and the first order inflaton perturbation S^tp are not 
independent. Once either <j>^ or S^'tp is known, the 
other may be computed from l|15|) . though it is simplest 
to solve for (f>^ and use (|15|) to compute S^'ip. One 
obtains a dynamical equation for the metric perturbation 
by applying (EJ to the sum of S^GR = K 2 S^T§ and 
S^G] = k 2 S^T- Einstein equations. The result is 



-2[n 



2 in' - n 



(16) 

Notice that the perturbed Klein-Gordon equation for the 
inflaton is not needed to close the system of equations. 
We discuss the slow roll solutions of (fTB|l in some detail 
since later on we will need to construct the Green func- 
tion for the operator in the left-hand-side of l|l(jf> . 
The conformal time slow roll parameters are 



6 = 1 



e-r) = 



Wo 



7V_ _ K^ifl 

n 2 ~ 2 n 2 ' 
- i. 



(17) 
(18) 



During a pure deSitter phase e = r\ — and the scale fac- 
tor evolves as a(r) = —l/(Hr) with H constant. During 
inflation, however, the Hubble scale evolves slowly as i|17|) 
so that for small e one has 



o(r) 



1 1 

HtI- 



and 



H = — = a = aH = — — 
a t(1 



(19) 



The dynamical equation for (|16[) can be rewritten 
in terms of the slow roll parameters as 



r (77-2e) + fc 5 



^ = (20) 



where we used l|19[) and dropped higher order terms in 
e, rj. Treating the slow roll parameters as constant, the 
equation l|2(J|) has an exact solution 

= (-r) 1 ^-' [ Cl (k)HP(-kT) + c 2 (k)Hi 2 \-kr) 

(21) ' 

where v = 1/2 + (3e — rj) to lowest order in slow roll 
parameters. 

It remains to fix the coefficients ci(fc), C2(k) in Q2ip. 
The variable in terms of which the action is canonically 
normalized is the Mukhanov variable Q 



where the inflaton fluctuation is solved for using the con- 
straint equation l|15l) . We fix Ci(fc), C2(k) by requiring 
that = e~ lkT /V2k on small scales —kr » 1 which 
corresponds to the usual Bunch-Davies vacuum choice. 
This leads to 



Cl W=^/^- P[ ^ + 1/2)] 



M p V 2 2k 



c 2 (fc) = 0. 



The solution for the metric perturbation, then, is 



f) = /^ A exp[ f( , + l/2)] 
fe w V 2 M„ 2fc 



X (-^/^-^(-fcT). 



(22) 



It is straightforward to construct the inflaton fluctuation 
using (TSJ. One may also compute the comoving cur- 
vature perturbation at first order, = TiVS 1 / '{atp' ) , 
and verify that the solution (|2~2"ll reproduces the usual in- 
flationary prediction for the power spectrum. One may 
also verify that in the limit e — > and rj — > the analy- 
sis of this subsection reproduces deSitter mode functions, 
which are discussed in appendix A. 



Note that during a slow roll phase we can treat e and 77 as 
constant (even though H is not exactly constant) since 
e', 77' are second order in slow roll parameters: 

e' = -2e(77 - 2e)H < eH, rj = 2e rfli < rjH. 

The above equality for 77' is only strictly correct if 
d 3 V/dip 3 = which is true in the case of interest. The 
statement that the slow roll parameters are approxi- 
mately constant is, however, quite general. 2 



C. Conditions for a slowly varying tachyon mass 

Now we come to an important point for this paper, 
that if rj is sufficiently small, then the tachyon is also a 
light field during some part of the observable period of 
inflation. The tachyon mass is given by m 2 = —Aw 2 + 
g 2 <^Q(/j). If we choose the arbitrary time t s in (|14|) to be 
when m 2 = 0, then g 2 ip 2 = Aw 2 , and 



ID' 



Of course this statement generalizes to cosmic time as well. 




-^j I -2r)\v 2 H(t-t s 



(23) 
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where N is the number of e- foldings of inflation occurring 
after the tachyonic instability begins. At some maximum 
value N — N*, inflation will end. If the inflaton rolls 
slowly enough, then the tachyon mass remains close to 
zero for a significant number of e-foldings. In general, we 
will have — iVj = N e — e-foldings of inflation before 
the spinodal time, followed by iV» e-foldings during the 
preheating phase. 

The approximation in (|23|l . that the tachyonic mass 
changes slowly enough for its time dependence to be 
approximated as linear, is true so long as |2rjiV| <C 1. 
This can be rephrased using the definition of rj in (|f 61 . 
and eliminating mn, using the COBE normalization of 
the inflationary power spectrum: V/{M^e) = 1507r 2 (2 x 
I0~ 5 ) 2 = 6 x 10~ 7 (see for example ref. [I]]). Using 
V — jXv 4 and eq. (|f 3(1 for e, the COBE normalization 
gives 



ml = 230 gX 



M 3 



(24) 



Then with N ~ 60, the requirement \2r)N\ < 1 becomes 



«= M v 
g < 10" 5 — ^ 

v 



(25) 



Interestingly, the bound i|25[l turns out to be a require- 
ment that must often be satisfied for different reasons, 
namely the experimental limit on the spectral index of 
the first-order inflaton fluctuations. In terms of the slow- 
roll parameters, the deviation of the spectral index from 
unity is given by 

n — 1 = 2r\ — 6e = 2r/ 

where we have used the fact that e <C rj in hybrid inflation 
(this is equivalent to the requirement that the energy den- 
sity which drives inflation is dominated by \v A /A). The 
experimental constraint on the spectral index is roughly 
\n — 1| ;$ 10" 1 . Writing n in terms of model parameters 
this translates into the constraint 
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5 x 10" 



(26) 



This is just five times weaker than the technical assump- 
tion (|2ii|) . 

In passing, we note that there is also a lower bound 
on g from the assumption that the false vacuum energy 
density is dominated by Au 4 /4 > rr?Ap\l2. Using g 2 ip 2 = 
Xv 2 and (|2"lfl . one finds 



g > 460 A —pr 

M 3 



(27) 



for the tachyon field fluctuation in Fourier space, 



6^ 



6H5 {1 



CTfc 



~n2 



112 

[9 <^o 



Xv 2 



6' 1 



CTfc = 0. 

(28) 

Once ipo < X x l 2 vj g the tachyon effective mass parameter 
d 2 V/da 2 becomes negative and the fluctuations S^a are 
amplified due to the spinodal instability. The efficient 
transfer of energy from the false vacuum energy Au 4 /4 to 
the fluctuations S^a is refered to as tachyonic preheating 
in the literature p^-p^. 

The initial studies of tachyonic preheating focused 
on the flat space dynamics of S^ak in the instanta- 
neous quench approximation. In this approximation the 
field S^'a is initially assumed to have zero mass and at 
t = a negative mass squared term — \m 2 \(8^ a) 2 / 2 is 
turned on. We briefly review these dynamics here fol- 
lowing closely p3|. Initially the tachyon has the usual 

Minkowski space mode functions e~ tkt+lk ' x /V2k 3 . Once 
the negative mass squared term is turned on the modes 
with k = \k\ < \m a \ grow exponentially with a dispersion 



(S^af) = -L /""" dkke 2t V^F^ (29 ) 



\m„ 



An 2 a 



^(e-(a-l) + l); 



a = 2 1 m a \t 



which produces a spectrum with an effective cutoff 
^max = \ m a\- The tachyonic growth persists until the 
dispersion saturates at the value 



(<5« 



1/2 



(30) 



at which point the curvature of the effective potential 
vanishes and the tachyonic growth is replaced by oscil- 
lations about the true vacuum. This process completes 
within a time 



1 ,n' 



2 mJ V X 



(31) 



which we call the spinodal time. At this point a large 
fraction of the vacuum energy Xv J A has been converted 
into gradient energy of the field 5« a so that the universe 

1/2 

is divided into domains with ({S^ <r) 2 ) = ±K of size 
I ~ |r7i CT | _1 and on average one still has (a) = so that 
o~o(t) — 0. These analytical argument are backed up by 
semi-classical lattice field theory simulations in j2^, • 
The discussion of the dynamics of tachyonic preheat- 
ing above apply strictly only in Minkowski space. The 
dynamics of tachyonic preheating including the dynamics 



D. Tachyonic Instability 

To quantify the evolution of the tachyonic instability 
at the end of inflation, we consider the equation of motion 



3 Even in deSitter space the Bunch-Davies vacuum choice will en- 
sure that this behaviour is respected on small scales k S> aH 
regardless of the tachyon mass during inflation. See appendix A 
for a review. 
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of the inflaton but neglecting the expansion of the uni- 
verse were considered in 24] . The dynamics of tachyonic 
preheating including both the dynamics of the inflaton 
and the expansion of the universe were considered in |35T ] 
wherein the authors reach conclusions identical to those 
discussed above. The authors of also find that the 
spinodal time is somewhat modified from 1|3I[I due to the 
background dynamics (see also (36L l37| ). 

Notice that we do not need to replace the average of 

1 /2 

the fluctuations ((<5^V) 2 ) with an effective homoge- 
neous background cr (r) because we work to second order 
in perturbation theory and the effect of these fluctua- 
tions enters into the calculation through the second or- 
der perturbed energy momentum tensor (S^'Tfy. When 
(<5( 2 )T^) becomes sufficiently large the backreaction will 
stop inflation. We take (|3U|) as our criterion for the end 
of inflation. We have checked numerically that this is 
a somewhat more stringent constraint than demanding 
that the energy density in the fluctuations S^a does not 
dominate over the false vacuum energy which drives in- 
flation Xv 4 /4. 

In the present work, we are interested in a situation 
which is different from the instantaneous quench, where 
the instability may turn on slowly compared to the Hub- 
ble expansion, rather than suddenly. We are approximat- 
ing the time dependence of the tachyon mass as being 
linear around the time when it vanishes, eq. (|23|) . so the 
mode equation can be written in the form 



prescription which is discussed in subsection IIII Bl and 
also in appendix A. 

Since (|32[l has no closed-form solution, we approxi- 
mate it in two regions. First, when ke~ 2N > cjJV|, 
we ignore the mass term and use the massless solutions, 
£k — a~ 3 / 2 -ff'y \{ke~ N ). We match this onto the solu- 



fc ■ - w "3/2^ 

tion in the region where k 2 e~ 2N < c\N\, where the term 
k e~ 2N is ignored in the equation of motion. The tran- 
sition between the two regions occurs at different times 
iVfc for different wavelengths, given implicitly by 



N k = In ■ 



k 



ln^Vfc 



(36) 



This is a multivalued function of x = k/y/c, because for 
x < (2e) -1 / 2 , the k term in the differential equation 
comes to dominate again for a short period around N — 
0, the moment when the tachyon is massless. To deal 
with this, we are going to assume that the solution is 
still well-approximated by the massive one during this 
short period. This amounts to replacing the multivalued 
function with the single- valued one shown in figure ^ 
We checked this approximation in conjunction with other 
approximations we will make for the mode functions, as 
described below eq. (|3"5|l . 



dN 2 dN k 



k 2 e- 2N - cN 5™a k = 0. 

(32) 

where TV = H(t — f s ), k = k/H and, incorporating the 
COBE normalization as in (|24(l . 



22000 gM p /v. 
From eq. (|25|l . c is limited to values 



c < 



M r 



(33) 



(34) 



for the validity of the approximation that the tachyon 
mass squared varies linearly with time. The quantum 
mechanical solution in terms of annihilation and creation 
operators a k , at has the usual form 4 



d 3 k 
(2tt) 3 / 2 



a k £ k (N) e 



■ikx 



h.c. 



(35) 



but the mode functions £fc will be complicated by the 
time-dependence of the tachyon mass. We normalize the 
mode functions t; k according to the usual Bunch-Davies 



N, 




single-valued approximation 



0.5 



FIG. 1: Exact solution and our approximation for the function 
N k in eq. l!^ . 



In the second region, with N > N k , the solutions are 
approximated by Airy functions, but it is more conve- 
nient to use the WKB approximation to obtain an ex- 
pression in terms of elementary functions. Ignoring over- 
all phases, in this way we obtain 



4 Our conventions for fourier transforms and mode functions are Cfe 
discussed in detail in appendix D. 



, -I/ 2 

2Hk 3 I I 



ike 



—N 



N <N k 



"(I + kl)- 1/4 , N>N k 



(37) 
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with 



and 



/2^( C |^l) 3/4 exp(£^ /2 ) 



z= 1 



-dV 



^fe = 1 
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c7V fc 



(38) 



(39) 



In the above expressions, we have for simplicity matched 
the amplitudes but not derivatives of the solutions at 
N = Nk- This will not affect the estimates we make be- 
low. We used to reexpress exponential dependence 
on Nk as power law dependence. Notice that exponent 
in i|38|) becomes purely imaginary when |ciV/. < — 1. We 
also replaced Iz] 1 / 4 — > (1 + l^l) 1 ^ 4 to correct the spuri- 
ous singularity at z — where the WKB approximation 
breaks down. We numerically verified that this gives a 
good approximation to the exact Airy function solutions. 

Moreover, we have checked the approximate solution 
by numerically integrating the mode equations, starting 
from the small- TV region k 2 e~ 2N ^> c\N\, where the mass- 
less solutions with known amplitude tell us the initial 
conditions, and integrating into the large- N region where 
the exponential growth due to the tachyonic instability 
becomes important. We did this for two orthogonal solu- 
tions to the mode equations, £1,2, whose behavior in the 
small- N region is 

& = {2Hk 3 )- 1 e~ N (kcos(ke- N ) ~sm(ke- N )) 

Ca = {2Hk 3 )- 1 e~ N (ksm(ke- N ) +cos(ke- N )) 

Evolving these initial conditions to large N, the enve- 
lope of these functions, whic h is also the modulus of the 
complex solutions, is £ = + £f • We compared this 
numerical solution to the modulus of our approxima- 
tion l|37|) for a large range of c and k values. In the large 
N- region, agrees with £ up to a numerical factor of 
order unity. This factor, the ratio of the actual solution 
to the approximation, is shown in figure (J2J). Because the 
exponential growth of the mode function is a very steep 
function of N, these small errors have an imperceptible 
effect on the exclusion plots we will present in section IVll 
Furthermore, we have checked which values of c and k ac- 
tually give constraints in the parameter space of the hy- 
brid inflation model below, and found that in the regions 
where c is large, k is exponentially small. Extrapolating 
the results of figure J5J) indicates that the error becomes 
quite small as k — > 0. Therefore our approximations for 
the mode functions are quite good for the purposes of 
this paper. 

With the above approximate solution, we are in a po- 
sition to recompute the dispersion of the tachyon fluctu- 
ations in the case of a more slowly varying tachyon mass, 
((S^a) 2 ) = (2tt) 3 J d 3 k\£ k \ 2 . Following the discussion 
of subsection 1111 Dl we set this equal to u 2 /4 at the end 
of inflation, N = N*: 



d 3 k 
(2tt) 3 



(40) 



1.4 
1.2 



+ 
_ 


1 1 1 1 


< • c = 0.005 

c = 0.01 

c = 0.1 

c=l 
c= 10 
c= 100 




1 1 1- 

.1 

I'll i 
: ]\ » 

! \ i"\ 








— 


\ ^ - 


- - " 


; f"'^ 1 ) 

1 




\ s" 




1 




s 

„ A 


x 


, >* 




_ - " \ 


/' \ 






" \ 

\ J 


'■ \, , 






1,1, 


i \ 


1 



-2 



-1 1 

log, n km 



FIG. 2: Ratio of the exact mode functions to the approxima- 
tion 13711 . at large times. 



which implicitly determines in terms of parameters of 
the hybrid inflation model, 



N*=N*(g,\,v/M p ) 



(41) 



The main contribution to the integral at N = comes 
from wave numbers for which the exponentially grow- 
ing solution in (|37|l applies. These modes satisfy k < 
k max = He N * \fcNl. We have numerically performed 
the integral for a wide range of values of c and N„. 
The result is displayed in figure [21 where contours of 
lnM 2 /Xv 2 are shown in the plane of and Inc. Re- 
call that c = 22000 g (M p /v), eq. f3"3"Jl. 




N* u 10 12 14 

FIG. 3: Contours of \nM^/Xv 2 in the plane of iV* and In c. 



N=N, 
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IV. SECOND ORDER FLUCTUATIONS IN THE 
LONG WAVELENGTH APPROXIMATION 

A. The Master Equation 

The authors of have derived a "master equation" 
for the second order potential (f)^ which can be written 
as 



2H( V 



b'W + [2-H 2 ( V -2e)-d k d k ]^ 
= J(t,x) (42) 



where the source terms are constructed entirely from 
first order quantities and can be split into inflaton and 
tachyon contributions 

J{t,x) = J a {T,x) + J v (t,x). 

Although we have explicitly inserted the slow roll param- 
eters, equation (|4*2)) is quite general and we have not yet 
assumed that 77, e are small (although recall that we have 
set <7o = 0). We have verified both the second order Ein- 
stein equations and the master equation presented in |32^ | 
and these results are discussed in appendix B. 

Because the equation (I42|) is linear we can split the 
solutions into three parts: the solution to the homo- 
geneous equation, the particular solution due to the infla- 
ton souce and the particular solution due to the tachyon 
source. The solution to the homogeneous equation will be 
proportional to (f)^ since the differential operator on the 
left hand side of (j42(l is identical to the operator which 
determines (j)^ , equation (|16[) . In cosmological perturba- 
tion theory the split between background quantities and 
fluctuations is unambigiuous, since background quanti- 
ties depend only on time while fluctuations depend on 
both position and time. However, the split between first 
order and second order fluctuations is arbitrary and the 
freedom to include in the solution for <f>^ a contribution 
which is proportional to <j>^ reflects this. We fix this 
ambiguity by including only the particular solutions for 
0( 2 ) which is due to the source, J(t,x). 

During the preheating phase the tachyon fluctuations 
are amplified by a factor v/H and so J a will come to 
dominate. This implies that the particular solution for 
^( 2 ) which is due to J a will come to dominate over the 
particular solution which is due to J v . Although our 
analysis will focus on this part of the solution, we will 
also consider the inflaton source in order to verify that 
our formalism reproduces previous results. 



equation, which is equivalent to assuming that the second 
order fluctuations are zero before the source is turned on. 
During a quasi-deSitter phase the master equation can be 
written as 



Or 



= (-t) 2 ^ 



-(r,-2e) + k 2 
Jfe(r) 



(4 



where it is assumed that e,r) <C 1, and the differential 
operator on the left hand side of the master equation 
is written in a manifestly self-adjoint form. In deriving 
B43JI we used Ijl9|l and treated the slow roll parameter as 
constant, which is consistent at first order in the slow roll 
expansion. The causal Green function for this operator 
is 



Gk (r,r') = |e(r-r')(rr') 1/2+ "- £ 



(44) 



X [J v (-kT)Y v {-kT') - J v {~kT')Y u {~kT)] 

where the order of the Bessel functions is v = 1/2 + 3e — r/ 
as in subsection IIII Bl 

The solution for the metric perturbation cfr- 2 ^ can then 
be written as 



dT'G k (T,T')(-T , f e - T > ) J k (T') 



-(l+e)/a t B 



where a% = a(U) is the scale factor at the some initial 
time, well before the tachyonic instability has set in. This 
solution is quite general; it applies during any slow roll 
phase, including during the tachyonic instability. 

There are several interesting limiting cases of <|44ll . In 
the long wavelength limit k — > the Green function re- 
duces to 

G fc (r,r') = e(r-r , )(l + 2r 7 -6e) 

_(_ r )l+2e /_ T A2(»j-2e) + (_ r m+2e/_ T \2(r)-2e) 

In the case of pure deSitter expansion e = r\ = 0, but for 
all k, the Green function reduces to 

G k (r, t') = 6(r - r')^ sin [k (r - r')] . 

Finally, the form of the Green function in the case of 
e = r\ = and k — > may be of some interest 



G k {T : T') = Q{T-T'){T-T'). 



(46) 



B. The Green Function 

In this subsection we construct the Green function for 
the master equation so that (f)^ may be determined in 
terms of first order quantities. As discussed previously 
we consider only the particular solution for (j) 1 ^ due to 
the source J and neglect the solution to the homogeneous 



C. The Gauge Invariant Curvature Perturbation 

The curvature perturbation is expanded to second or- 
der as 

c = c (1) + ^c (2) 



9 



In hybrid inflation, where oo = 0, the first order contri- 
bution comes entirely from the inflaton sector 



H 



5W, 
P'o 



where p = -(T °) (0 ), <5 (1) p = -<5 (1) T ° are the unper- 
turbed and first order stress tensor respectively. It is 
also conventional to define the comoving curvature per- 
turbation at first order 



<P'o 



which, on large scales, is related to £W as 



rc« + c (1) = o. 



The definition of the first order curvature perturbation 
is generally agreed upon in the literature (up to the sign 
of C^ 1 "*)- At second order, however, there are several def- 
initions of the curvature perturbation in the literature 
(see 0| for a comprehensive discussion). The definition 
4(1 and generalizes the definition of 
ill (valid on large scales) 



we adopt follows 
Malik and Wands 



C (2) = V^ (2) + ^ {2) P 2 7 %8^p8^ P ' (47) 
Po (Po) 



Po V 1 



J\2 



Po 



H' 



- 2 



to multiple scalar fields. The definition of the curvature 
perturbation adopted in j^] generalizes @ to two scalar 
fields. However, it has been shown that this definition 
applies only during inflation 0, 0] since it is conserved 
only in the slow roll limit. 

Because Co = the inflationary trajectory is straight 
the case of hybrid inflation and £W is conserved on large 
scales l4 ^ -|45 | . However, this is not the case at second 
order and one expects that will be amplified due 
to the tachyonic instability. This amplification of Q 2 ^ 
will continue until N = N* at which point the backreac- 
tion sets in and stops inflation. For N > iV* the inflaton 
is no longer dynamical since all of its energy has been 
converted into tachyon fluctuations. Thus for N > 
the large scale curvature perturbation is conserved at all 
orders in perturbation theory since only one field (the 
tachyon) is dynamical and there are no non-adiabatic 
pressures. 

In jijj the second order large scale curvature pertur- 
bation is written in terms of the first and second order 



Sasaki-Mukhanov variables as 
1 1 



£(2) = 



+ 



3-e(^ 
1 1 

1 1 

3 - e W 

4(3 + e - n) 



v'oQT 



Q 



'(i) 



y(i) 



Q 

3- 2e 
3- e 



2 2 

a m„ 



i 2 2 

+ a m v 



(«<«y 



H 
— < 



,(i) 



+ (-10 + 26 + 2??) (— Q« 



H 



(48) 



where the first order Sasaki-Mukhanov variables 5 are 
QW = (49) 
Qi 1} = <5 (1 V (50) 



and the second order Sasaki-Mukhanov variable is 



Q< 2) 



<^(2) 

H 



+ (2 + 2e-r ? )^(^)) 



(51) 



+ 2^y« + A^v. 

In writing (|48I51|) we have restricted ourselves to hybrid 
inflation, made use of the background equations and in- 
serted e,r] (but not assuming slow roll). 

The last four lines of l|48|l are relatively simple to eval- 
uate in the large scale and slow roll limit. On the other 
hand, the first line of l|48|) is somewhat more compli- 
cated and its evaluation requires the second order Ein- 
stein equations. Thus we focus here on the following 
contribution to <^ 2 ) 



c 



(2) 



1 



3-e(^) 2 



M 2) 



(2) 

From the definition of Qlp ' it is clear that this contains a 
contribution of the form (p' 5^ 2 'ip' + a 2 dV/dip5 l - 2 ' > ip which 
also appears in the second order (0,0) Einstein equa- 
tion (|B-1J| . We therefore use l)B-lj) to eliminate S^ 2 'ip 
from 15211 . We also eliminate tp^ in favour of 4>^ using 
(|B-lip . The result is that J53J) takes the form 

1 d k d k 4>W 



dip 



,(2) 



(52) 



(2) 



3 - 



eH \e 



1,(2) 



1 

7R 
1 



Ti 1 



3-e 
1 1 




m 2 



3-eeH 2 



3 



3 - e eH 2 
T' 



a 



dV, 



T 4 (53) 



5 Notice that is related to the variable yW discussed in sub- 
section HSU by Vt 1 ) = aQ^K 
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where the quantities Ti,7 are constructed entirely from 
first order fluctuations and are defined explicitly in ap- 
pendix B in equations {El [ED EH USE ESI EH- 
The quantity T4 is also contructed from first order fluc- 
tuations and is defined as the last two lines of (l5ll) . That 



The quantity 7 CT can be written in the form (see equation 
(66) of [13|, or equuivalently 



la 



-K A" 



is, 



T 4 = 



(2 + 2e-r?)^ 



(54) 



39, (d k d k 5 {1) ad l 5 {1) <7 



From l|53|l we see that the dependence of on the 
second order fluctuations comes from the combination on 
the first line of l|53[) which can be computed in terms of 
the source using the Green function (|44H . 

We split £( 2 ) into contributions coming from the infla- 
ton and the tachyon as 

c (2) = c< 2) + d 2) 

and study each piece separately. This splitting is differ- 
ent from the one discussed in [4(J], where the curvature 
perturbation is defined for each fluid in such a way that 
the total curvature perturbation is a weighted sum of the 
individual contributions. Instead, we simply divide £( 2 ) 
into terms which depend respectively on the tachyon and 
inflaton fluctuations, 5^ a and ^ , which is only 

possible because cto = 0. 

(2) 

The tachyon part of the curvature perturbation Q 
gets contributions from the first and second line of 
()48p. both explicitly through and implicitly through 

S^(p, . The inflaton part of the curvature perturba- 

.(2) 



Notice that the terms in the first line of the l|55|l are lo- 
cal, the terms in the second and third line are non-local 
(containing an inverse laplacian A -1 ) and the fourth line 
contains terms which are both local and doubly non-local 
(containing A~ 2 ). The Fourier transforms of the source 
terms are computed in appendix D wherein we also dis- 
cuss our conventions for the inverse laplacian operators. 

In the following, we will need the Fourier transform of 
terms like A 7 CT , 



d 3 k' 
fc 4 J (2tt) 3 / 2 



k' 2 k-(k-k')5^& k ,6^& k _ k , 



K 

2k 2 



d 3 k' 
(2tt) 3 / 2 



fc'-(fc-fc')(5 (1 W(5 (1) CTfc-fc< 



This expression is operator-valued and can be written in 
terms of annihilation/creation operators and mode func- 
tions as 8^a k = a k ^ k {t) + a_ fe £_fc(i) (see appendix D 
for more details). In the Fourier transformed expression 



tion QT contains contributions from the last three lines for A la, the scale dependence of the mode functions 

S^'ak is integrated over so that A _1 7 CT gets contribu- 
tions from the tachyon fluctuations on all scales. The 
large scale limit of terms like A~ 1 j (7 is not transparent 
and therefore we do not neglect any terms in the tachyon 
source which contain inverse Laplacians. 

>(2) 



of (|48|l as well as from the first line of |@SJ, both implicitly 
through S^(p, and explicitly through the definition 
of Q {2) . 

The inflaton part of the curvature perturvation, 



(2) 

'■P 

coincides with the C^of single field inflation and has 

been derived previously 0, El 03 • We have considered 

(2) 

the construction of Q, using our formalism and these 
results are presented in appendix C. 



D. The Tachyon Curvature Perturbation 

We now consider the tachyon contribution to the sec- 

(2) 

ond order curvature perturbation Q- . It is sourced by 
J a which, in position space, takes the form (see equation 

r(T,x) = aVm 2 (>V) 2 -2 K 2 (>V) 2 

+ 2k 2 H(1 + 77 - e)A- 1 d l (S^a'd^^o) 

+ 4k 2 A- 1 9 t 5 4 (S^a'd^a) 

- H(l + 2e-2 V )A- 1 1 'a + &- 1 l"- (55) 



We now compute Q . The tachyon contribution to 
^ 2 ) comes about entirely through the implicit depen- 
dence of <5^ 2 V and ip^ on 5^a in the first line of (gSJ 
and the explicit dependence on the second line of 14811 . 
Our focus is on the leading order contribution to in 
the slow roll and large scale limit. If we work only to 
leading order in the slow roll parameters it is sufficient 
to use the Green function (|44|) in the limit e = r\ = and 
keep only the terms in the tachyon source (|55|l which are 

not slow-roll suppressed. However, to consistently com- 
(2) 

pute £<r we must keep the next-to-leading order terms 
in the small (k/TL) 2 expansion of the Green function. To 
see this notice that powers of k 2 cancel inverse Laplacians 
in the source: 



k Jl — 



~la,k 



ni 



cr,k 



where • • • denotes gradient terms which are small on large 
scales. In appendix B it is shown that 7^ can be written 



11 



2 -a 2 ml (W* 



on large scales as (|B-14|) 



3k 2 
~2 



Thus k 2 J£/1i. 2 contains terms which are of the same 
form as those which appear in J£ and hence these 

(2) 

terms must be included to consistently study Q on large 
scales. 

The curvature perturbation (^ 2 ) depends on 0( 2 ' 
through the combination —(fy'^/eTt — <j)^/e + 
(9 fc <9 fc 0( 2 ) /3eTi 2 at leading order in slow roll param- 
eters Using the Green function for e = r) = 
expanded up to order k 2 we find that 



1 '° 

e 



4 2) 



e 3eH 2 
dr'9(r-r')^(r') 



1 



' 2 I 
-T T 



(56) 



Integrating l|56[l by parts 6 and plugging the result into 
(|48ll and l|53l) we find that on large scales 



'(2) 



^') 2 
H(r)3 



(VV) -a 2 m 2 



(57) 



where the tachyon fluctuations S^a are functions of the 
integration variable r'. The corrections to ()57f) are either 
subleading in the slow roll expansion or are total gradi- 
ents which can be neglected on large scales. Equation 
(|57|l is the main result of this section. The interested 
reader may find a more detailed discussion of our calcu- 

(2) 

lation of Q in appendix E. 

It is interesting to contrast the simplicity of this result 
with other previous computations of ostensibly the same 
quantity |22j . Part of the difference is due to our different 
definition of the second order curvature perturbation, but 
one must also take considerable care in keeping all terms 
which are of the same order in slow roll parameters and 



powers of k 2 in order to obtain the delicate cancellations 
that collapsed the sum of many terms down to this com- 
pact form. Especially notable is the absence of nonlocal 
terms (containing inverse Laplacians) in the final result, 
which were quite prevalent in the intermediate steps. It 
is an important consistency check that such terms dis- 
appear in the end, since they do not respect causality. 
For example, a term of the form A~ i diS^ 1 'ad i S^ 1 'a gives 
rise to an acausal response to a source which is local- 
ized in time and space. Consider a source of the form 
S^ 1 ' a ~ e~ x / a which is turned on at some instant in 
time. Then A _1 9j#W ad l S^ a ~ l/|x| instantaneously 
at large distances, instead of being exponentially small. 
This would be a clear violation of causality and is phys- 
ically inadmissible. 

To clarify, we have shown that at leading order in slow 

(2) 

roll parameters the nonlocal contributions to Q cancel. 
We have not checked that such terms cancel at higher 
order in the slow roll expansion, though we believe that 
they do. This cancellation was not observed in previous 
studies because the subleading (in k 2 ) corrections to the 
Green function were not included and thus the large scale 
expansion was inconsistent. 

A comment is in order concerning the long wavelength 
approximation which we have used in deriving (|57l) . In 
writing the expression for k 2 J£ we have dropped terms 
which are total gradients even though such terms will be 

(2) 

integrated over time in computing Q , due to the time 
integral in l|45|) . Strictly speaking, one should only ap- 
ply the large scale limit kr -C 1 after the time integral 
has been performed. The reason for this is that the time 
integral extends from when the modes are well within 
the horizon, where they oscillate, to when the modes are 
outside the horizon, including horizon crossing. Mode- 
mode coupling near the epoch of horizon crossing can 
contribute momentum-dependent terms which are not 
suppressed on large scales [jj, |5fJ. However in our ap- 
plication, we are justified in dropping the contributions 
which are generated near horizon crossing since they are 
exponentially suppressed compared those which are pro- 
duced during the preheating phase. 



E. Time-integrated tachyon perturbation 

We can make our main result (|57|l more explicit by sub- 
stituting in the solutions I|35I37|I . The result is simplified 
by taking the Fourier transform of C*- 2 ) and evaluating it 
at vanishing external wave number, and at the final time 
corresponding to the end of inflation, N = TV* : 



This integration by parts will give rise to terms which are 
evaluated between the initial and final time of the form 
[• • • ]t'=— 1/ojH' Since we are interested in the preheating phase 
during which the fluctuations are amplified exponentially we can 
safely drop the contribution at r' = —1/a.iH relative to the con- 
tribution at t' = r. 



d 3 p 



( a p bp a P b- p + perms) 



K 

e J (2tt) 3 / 2 

cN, 

\ / f{c,N,N*)dN (58) 



<N p ,Ni) 
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where iVj denotes the value of N at the beginning of 
inflation, /(c, AT, TV*) is given by 

/( C ,iV,iV <f )=e- 3JV +^ 3/2 (l + Hr 1 / 2 x 



4 V 



3(N-N,) 



1-e 



ciVe 3(^-iV, 



/i-1- 



2c sign(z) 



27(1 + |*| 



(59) 



where £ = (1 + |c/V), and "perms" in l|58[l indicates the 
three other combinations of a p b p and a* p 

For illustration, we show the behavior of the function 
f(c,N,N*) for sample parameter values N* = 22.5 and 
lnc = 1. Figure 0] plots (sign(/)ln(l + |/|) as a func- 
tion of N. The function is exponentially peaked at the 
initial value N = Ni, and at the final value N = N*. 
Moreover it always becomes negative at N* because the 
negative mass squared term in <|59|) comes to dominate. 
Although it is not obvious in the figure, the negative value 
is orders of magnitude larger than the positive maximum 
just preceding it, so the negative extremum dominates 
in the integral in (|59|l . Whether the extremum at AT* or 
Ni dominates overall depends on how |3ATj| compares to 
9/(2c)z 3/2 evaluated at AT*. 

Because of the exponential growth of / at its extrema, 
it is a very good approximation to the integral to write 
/ = e 9 and expand g = g m +g' m (N — N m ) in the vicinity of 
the maximum value, whether it is at Ni or A*. Since the 
integral is so strongly peaked near the extremum, there 
is an exponentially small error in extending the range of 
integration to the half-line. In this way one obtains 



dNf 



\!)n 



(60) 



We will use this approximation below to numerically eval- 
uate the integral. 



V. BISPECTRUM AND SPECTRUM OF 
SECOND ORDER METRIC PERTURBATION 



200 




FIG. 4: (sign(/) ln(l + |/|) versus N, showing behavior of the 
the function / defined in tSUI . for lne= 1, JV* = 22.5 



Carrying out the contractions of pairs of fields which 
contribute to the connected part of the bispectrum, one 
finds eight terms, which in the limit of vanishing external 
wave-numbers, are all equal. The result is 



B = 



d 3 p 

(27T) 5 



\b P \ 



max(JV p ,iV i ) 



dN f(c, N, N* 



(63) 



The integrand of the p-integral is exponentially strongly 
peaked, either near \a(p/^/cH) = N{ if / has its global 
maximum near A" = N, or else near hi(j>/^/cH) = if 
/ is dominated by its behavior near AT*. The logarithm 
of the integrand for a typical case is shown in figure 
If the cusp-like peak is dominating, we can use the same 
approximation as in l|60l) to evaluate the p integral, on 
both sides of the maximum. If the other local maximum 
at larger values of p is the global maximum, as sometimes 
happens, then we should treat the integral as a gaussian, 
since the derivative of the integrand vanishes at the max- 
imum. In this case we similarly write the integrand of the 
p integral as / = e 9 and expand g = g rn + ^g'^ip - p m ), 
where g'^ < 0; then ll6TJ|) is replaced by 



Here we calculate the leading contribution to the three- 
point function (bispectrum) of the second-order curva- 
ture perturbation due to the tachyon, 



a(2) A2) 



A2) 



(2tt)- 3 / 2 S^ih +k 2 + fc 3 ) B(k u k2,k 3 ) 

(61) 

It is understood that only the connected part of the cor- 
relation function is computed, which is equivalent to sub- 
tracting the expectation value of ^ 2 -* from the quantum 
operator. This three-point function is straightforward to 
compute, using the free-field two-point functions 



(62) 



dNf 



' 2tt 



(64) 



We have carried out the evaluations numerically over a 
range of values of lnc and A^. In figure © we plot 
contours of In \B\, where B is defined by 



B = 



' " y f -X B 



2tt 2 



cc 



(65) 



Although we evaluated B at vanishing external wave 
numbers, fej = 0, the result is equivalent to evaluating 
B at very small wave numbers closer to horizon crossing, 
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FIG. 5: Log of the integrand of the p integral in eq. 16311 . for 
a case where the maximum occurs near \n(p/^fcH) = Ni . 
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FIG. 6: Contours of In \B\, defined in I65II . in the plane of N* 
and Inc. Ni — —30 in this example. 



hi = He *. This is not immediately obvious, but one can 
see that at very small p, the p integral in (|63[l behaves like 
J d 3 p/p 9 , which in the presence of nonvanishing external 
wave numbers would be J d 3 pp~ 3 (p — k)~ 6 ~ fc -6 . (This 
claim is explicitly verified in appendix F.) However, l|63|) 
does not diverge as k — > 0; the p integral is cut off near 
He * because of the lower limit max(A/j, N p ) in the inte- 
gral over N. If the lower limit were simply N p instead of 
max(ATj, N p ), the value of the N integral near its lower 
limit would go like e 3Np ~ 1/p 3 instead of e Ni . Thus the 
factor of e eN * in |JBSJ| is a reflection of the fact that the 
bispectrum (at least the part of it which is due to the 
low-p part of the integral) is scale- invariant, B ~ 1/fc 6 , 
for k > He Ni , but the infrared divergence is cut off for 
k < He Ni , and the spectrum remains flat as k — > 0. 
A consequence of the above discussion is that the bis- 



pectrum is scale-invariant only in the part of parameter 
space where the low-iV part of the N integral dominates, 
namely the lower left-hand region of fig. 10- The upper 
right-hand region of fig. JSJ gets its dominant contribu- 
tions from the large-p part of the p integral, which is not 
infrared sensitive. The bispectrum really is flat as a func- 
tion of k, for small k, in this region. Furthermore, B is 
negative in this region, whereas it is positive in the other 
region. This means there is a very narrow, fine-tuned 
curve along which B vanishes, seen as the borderline be- 
tween the two different behaviors. This is a region which 
will always be allowed as long as no nongaussianity is ob- 
served in the CMB, but since it is a set of measure zero 
in the parameter space, we will not take it into account 
in deriving limits below. 

The boundary between the two regions in fig. JSJ is 
described analytically by the relation 



m\ = Yc * 



3/2 



(66) 



with = 1 + |ciV*. This where the two terms in the 
exponent of / in (|59|) . which determinant the dominant 
behavior of /, just balance. In the region to the right, 
<i „3/2 ^ 3|jVj|, which is also where B is indepen- 
B is also independent of Ni. In the region to 
the left, B depends principally on Ni through the factor 
e -6Ni m (|65|^ However there is some mild extra depen- 
dence on Ni, as shown in figure Q. 



where w-z 
dent of k. 
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FIG. 7: Dependence of In \B\ on N f for AT* = 25 and lnc = -2. 
The smooth curve is the numerical fit In \B\ = — 8.3— 0.13iVj — 
0.0007JV, 2 . 

To make contact with experimental constraints, we 
want to compare the predicted bispectrum with that of 
single-field inflation, where nongaussianity is convention- 
ally expressed via a nonlinearity parameter Jnl, defined 
through 



B{h 



6 



fNL(P4>(ki)P <t> (k 2 ) + permutations) (67) 
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where P</,(k) is the usual inflationary power spectrum, 
P (fc) 1 / 2 - 1CT 5 (27r)/fc 3 / 2 . If we assume that all k,. 
the same, ki = k, then 



are 



IN 



(27r)- 4 10 20 P(fc)fc 6 



(68) 



We can thus convert our predicted bispectrum into an 
effective f^Li f° r which the present experimental con- 
straint is roughly |/jvl| < 100. To get the strongest 
constraint, k should be evaluated at scales near hori- 
zon crossing in the scale-invariant region of B, and at 
the smallest relevant scales in the region where B is fc- 
independent. Since B(k) is evaluated at the end of in- 
flation, the horizon-crossing scale is k ~ He~ *. On 
the other hand, current experimental constraints on non- 
gaussianity involve temperature multipoles of the CMB 
going up to I = 265 0| , which represents scales roughly 
e 5 times smaller than horizon-crossing. So we should use 
k = e e+5 i? in the region where B is constant at low 
k. It turns out that this is the region of parameter space 
where the experimental limit is saturated, so it is the 
relevant case. 

Using in (JHSJ and recalling that N e = |JV,-| + 
we find that 



Inl — 



-e 175 - 6N *c- 3 B 



(69) 



where we used H 2 = F/(3M|) as well as the COBE 
normalization V/(eMp) = 6 x 10~ 7 . The contours of 

B can thereby be converted into contours of /jvx ; an d 
demanding that \/nl\ < 100 gives a constraint in the 
parameter space c, , Ni . The limiting curves are shown 
for different values of iVj in figure |H1 These curves are 
quite insensitive to the actual value assumed for /jvx, 
because the effect turns on exponentially fast. The curves 
for /jvl = 1 are visually hard to distinguish from those 
shown (for Jnl = 100). 
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FIG. 8: Boundaries of the regions excluded by the nongaus- 
sianity constraint |/jvl| < 100, in the plane of lnc and iV*, 
for different values of Ni. 



A. Constraint from tachyon contribution to 
spectrum 

The analogous calculation can also be carried out for 
the tachyon contribution to the two-point function of the 
curvature perturbation. We show that it gives a stronger 
constraint than does nongaussianity. 

Closely following the preceding calculation, it is 
straightforward to show that, in the limit of vanishing 
external wave numbers, 



/•(2) >(2) 



2 K - 

i 



d 3 p 



= 5(ki + fc 2 ) x 



dNf(c,N,N*) 



^r p] j Np 

In analogy to ijMK]!. we define 7 

b - 4^ 2 c 3/2 IT ' 



(70) 



(71) 



and we display contours of S in fig. OH Analogously to the 
bispectrum, S is approximately scale-invariant (~ 1/fcf) 
in the left-hand region, but independent of k in the right- 
hand region. 




10 N 1 * 5 



FIG. 9: Contours of ln|S| in the plane of N, and lnc, at 
N = -30. 

The COBE normalization of the power spectrum im- 
plies that 



A2) A2) 



2tt 2 



< z TT 8{k l + k 2 )V c (72) 

con K 



7 The powers of H and c can be understood as follows: b p 
// _1</2 c~ 3 / 4 , whereas p ~ sfcH. 
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with -P c 1/2 = 10~ 5 . Combining this with J70|) and (jTTjl . 
we get an upper bound on the quantity 

where the strongest constraint is obtained by taking the 
largest k values which are measured by the CMB. We will 
conservatively take this to be k = e e He~ Nc . Notice that 
the tachyon fluctuations have a spectral index of n = 4 
(defined by k 3 S ~ fc™ -1 ) which is consistent with the 
Traschen integral constraints [H3 |. Following the same 
steps as we did for the bispectrum, this can be rewritten 
as 




0.5 1 1.5 2 2.5 

In c 

FIG. 10: Boundaries of the regions excluded by the spectral 
constraint < 1, in the plane of lnc and N*, for different 
values of JVj. The corresponding constraints from }nl (fig. El 
are shown as dotted lines. 

One might consider being more conservative and im- 
posing, say < 0.01 rather than < 1, as we have 
done. Our exclusion plots are actually quite insensitive 
to the value assumed for as was the case with /jvl, 
since the effect turns on exponentially fast. 



VI. CONSTRAINTS ON HYBRID INFLATION 
MODEL PARAMETERS 

In the previous section we computed the effective non- 
linearity parameter Jnl which characterizes the size 
of nongaussian fluctuations produced by tachyonic pre- 
heating in the hybrid inflation model under consider- 
ation. Experimentally, it is currently constrained as 
\Jnl\ ^ 100 [5l]], which is expected to improve to the 
level of \/nl\ < 5 from future experiments, like PLANCK 



[1113. Similarly, we defined a parameter _/x which char- 
acterizes the size of the second-order tachyonic contribu- 
tion to the spectrum relative to the experimental value, 
with the constraint < 1- These constraints were 
most easily determined in terms of derived parameters c 
and iV» , but here we want to recast them in terms of the 
fundamental parameters of the hybrid inflation model, 
the dimensionless couplings g and A, and the VEV of the 
tachyon field, expressed in the dimensionless combination 
v/Mp. We do not treat mt as an independent parameter, 
since we used the COBE normalization to eliminate it in 

Jill. 

The only obstacle to working directly in the model 
parameter space is the implicit dependence of on 
(A, v). We have numerically inverted the relation 
depicted in fig. © to determine N*(\,g,v), It then 
becomes straightforward to scan the model parameter 
space, recomputing Jl and Jnl- We display the con- 
straints in the plane of log 10 g and log 10 \ f° r a range of 
values of \og w {v/M p ). In performing this scan, one must 
keep track of whether various assumptions are respected. 
These include eqs. and (|T7Jl . respectively correspond- 
ing to the assumption that m 2 varies linearly with the 
number of e-foldings, and that the vacuum energy den- 
sity is dominated by Aw 4 /4 during inflation. In addition 
to these, we implement two others: we assume that (1) 
the reheat temperature exceeds 100 GeV, so that baryo- 
genesis can occur at least during the electroweak phase 
transition; and (2) we limit ourselves to considering val- 
ues c > 10~ 4 , which entails a lower limit on the coupling 
g depending on v, eq. Q4U[). This is a numerical limita- 
tion, due to the difficulty of evaluating the integral i|40|) 
for small c. However, analytic calculations can be done 
in this region, which we have carried out in Appendix F, 
as we describe below. 

In the region where the condition (|25|l . gv/M p < 10 -5 , 
is not satisfied our analysis breaks down because the 
tachyon mass m 2 varies exponentially rather than lin- 
early with the number of e-foldings. However, we do not 
expect that any significant nongaussianity will be gener- 
ated in the region where c is too large because in this case 
we will have had m^. > H 2 for a significant number of e- 
foldings before the instability sets in and the fluctations 
of the tachyon will have a large exponential suppression. 
However, to be conservative in our analysis we consider 
the region where 123[) is not satisfied as allowed by the 
experimental constraints on /l , $nl ■ This is a rather 
small sliver of parameter space, since the spectral index 
contstraint (|26l) requires that gv/M p < 5 x 10~ 5 . 

Once a set of values for A, g and v/M p are chosen and 
iV* has been calculated, one must still determine Ni in 
order to calculate the parameters /l and Jnl- We do so 
by first computing the total number of e-foldings using 
the standard result 

/10 16 GeV\ 1 ( V^ 4 \ 
^ = 62-ln(^^)--ln(^j (75) 

where V = Xv i /4 is the energy density during inflation, 
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and p r .h. is the energy density at reheating. In the fol- 
lowing, we will ignore the gravitino bound p~J^ <, 10 10 
GeV and assume instant reheating, p r .h. = V. The value 
of Ni then follows from Ni — N* — N e . We have checked 
that incorporation of the gravitino bound does not create 
a noticeable change in the excluded regions. 

The resulting constraints are illustrated for the case 
where v/M p = 10 -3 in figure ITT1 The unshaded region is 
excluded by our calculation of /l and /nl- The shaded 
region on the right falls outside our approximation H25|) 
or the spectral index constraint (|26() , that on top fails to 
satisfy i|27|) . and the region on the bottom has too low a 
reheat temperature. The region on the left has c < 10~ 4 
and we were not able to treat it using the same numerical 
methods, but we analyzed this region analytically in Ap- 
pendix F. The tachyon mass is small throughout inflation 
in this region, and gives rise to a scale-invariant spec- 
trum of nongaussian fluctuations, unlike the excluded re- 
gion shown in figure 1111 which corresponds to a scale- 
noninvariant contribution, with spectral index n = 4. 
Nevertheless, the magnitude of nongaussianity in this re- 
gion labeled, "also excluded," exceeds the experimental 
bounds. 
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FIG. 11: Unshaded region in the plane of log 10 g and log 10 A 
is excluded by the tachyon contribution to the second order 
curvature perturbation, for v/M p = 10~ 3 . Leftmost region is 
also excluded by a separate analysis (Appendix F). 

As was indicated in the previous section, we get a 
stronger constraint from the spectrum, whose distortions 
are parametrized by Jl l|73l) . than from nongaussian- 
ity, parameterized by /jvl ■ This is shown in figure 1121 
which plots the contours for the constraints /l < 1 and 
/jvl < 100 again in the plane of log 10 g and log 10 A, and 
for v/M p = 1CT 3 . 

The size and position of the excluded region is rather 
sensitive to the value of the VEV, v. We have varied 
log 10 v/M p between —1 and —9 to show how the con- 
straints change with v in figure 1131 One sees that the 
width of the excluded region for log 10 g grows linearly 




FIG. 12: Closeup of the borderline region in figure tTTI showing 
that the spectral distortion (fz, < 1, dashed line) provides the 
stronger constraint, relative to nongaussianity (/ml < 100, 
solid line). 

with log 10 Mp/ V i so the constraints are strongest at the 
smallest values of v. 




FIG. 13: Regions excluded by spectral distortion for 
log 10 f/Mp = — 1, — 3, — 5, — 7, — 9, as labeled on the figure. 
Note that for each v, the regions at smaller g are also ex- 
cluded, but due to nongaussianity (Appendix F) . 



VII. IMPLICATIONS FOR 
BRANE-ANTIBRANE INFLATION 

We now consider what implications our analysis has 
for the case of brane-antibrane inflation j5f|-[6(|, which 
is a string theoretic realization of hybrid inflation. In 
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this scenario inflation is driven by the interactions be- 
tween a brane and antibrane which are parallel to the 
three large dimensions and are initially separated in the 
extra dimensions of string theory. The brane-antibranc 
pair experiences an attractive force and move towards 
each other with inflation ending when the branes collide 
and annihilate. The role of the inflaton tp is played by 
the interbrane separation while the field a is the lightest 
streched string mode between the brane and antibrane, 
which develops a tachyonic mass when the separation be- 
tween brane and antibrane becomes of order the string 
scale. 

The potential of brane-antibrane inflation is not am- 
menable to the form for several reasons: 



• The inflaton potential typically has the Coulombic 
form V ~ A/(p n , rather than the simple form V ~ 

2 2 

• The tachyonic field is complex, rather than real. 



• The tachyon potential does not have its minimum 
at a finite value of the a field. 

• The tachyon field is described by a Dirac-Born- 
Infeld (DBI) action rather than a simple Klein- 
Gordon action as we have assumed in this paper. 

Here we will assume that these differences do not sig- 
nificantly alter the analysis and leave a more detailed 
investigation for the future. Specifically, this amounts 
to assuming that the behavior of the two models will 
be roughly the same if we match the quadratic parts of 
their actions, as well as the cross-coupling between the 
inflaton and the tachyon. Since the timescale for tachyon 
condensation is given by y— TOct , one might hope that 
the results are rather insensitive to the details of how the 
tachyonic instability saturates at large field values. 

In string theory the open string tachyon T between 
a Dp-brane and antibrane, separated by a distance y is 
described by the action 61] 




»Stac — ■ . 

C = V(T, y)yjl + Ms 2 g^d»Td v T^ 
where the small \T\ expansion of the potential is 

1 f(M s y) 2 r 



V(T,y) = r p 



2 V (2tt) 2 
G(|T| 4 ) + .. 



(76) 

(77) 
(78) 



In the above we have neglected the gauge fields and trans- 
verse scalars on the branes. The coefficient of the |T| 4 
term in the potential is order unity and its exact nu- 
merical value is unimportant for our estimates since the 
potential V(T, y = 0) is minimized only at infinite T. A 
potential of the form 



is typical. In the above M s is the fundamental string 
mass and t p is the Dp-brane tension. In order to obtain 
order-of-magnitude estimates for the model parameters 
A, g, v of © it is sufficient to consider only the tachyon 
part of the action l|76(l and neglect the inflaton poten- 
tial. We restrict ourselves to inflation models driven by 
D3-branes since inflation driven by higher dimensional 
branes have problems with overclosure of the universe by 
defect formation [fj^] . The D3- brane tension is given by 



T3 



(27T)3 ffs 

where g s is the string coupling. 



(79) 



A. The KKLMMT Model 

Brane inflation models in torodial compactifications 
are not realistic because the moduli which describe the 
size and shape of the extra dimensions were assumed 
fixed without specifying any mechanism for their sta- 
bilization. Recently there has been significant progress 
|58ll59| in reconciling brane inflation with modulus stabi- 
lization using warped geometries with background fluxes 
for type IIB vacua 63]. For our purposes the most im- 
portant feature of these compactifications is the presence 
of strongly-warped throats within the extra dimensions. 
The compactification geometry within the throat is well 
approximated by 



(80) 



ds 2 = a{yfg liv dx >l dx v + dy 2 + y 2 dttl 

where y is the distance along the throat, a(y) 
is the throat's warp factor and dfl 2 is the metric on the 
base space of the corresponding conifold singularity of 
the underlying Calabi-Yau space |64|. In the subsequent 
analysis we ignore the base space and treat the geometry 
as AdSs. 

A mobile D3-brane is placed near the UV end of the 
throat y = while an anti-D3-brane remains fixed at 
some location near the IR end of the throat y — yi > 0. 
The warp factor aj = a(yi) can be easily made much 
less than unity by a suitable choice of background fluxes 
and this large warping flattens the inflaton potential. 
This large warping also suppresses the tachyon poten- 
tial (effectively reducing the brane tension to afr^,). At 
quadratic order, the tachyon action takes the form 



8 = - 



x v-g 



atT3 



( M 2 y 2 
2 \(2n) 2 



<4n 
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2M? 



\8T\ 2 



(81) 



V(T,y = 0) 



\T\ 2 /a 2 



The effective values of the couplings can be found by 
rewriting the action l|81(l in terms of the canonically nor- 
malized fields a = a^pr^I '/M s and <p = ^Jr^y (see equa- 
tions 3.6, 3.10 or C.l in (53), and then matching to the 
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hybrid inflation potential ©. This gives the correspon- 
dence 

0; M. 



r 3/2 



/9s 



A = Y 9s 
g = a i y / 27rg s 



Before trying to make a general investigation of the pa- 
rameter space, we will examine the fiducial values which 
were considered natural in (58|: «i = 2.5 x 1CT 4 , g s = 0.1, 
t 3 /M 4 = 10~ 3 . The COBE normalization also fixes the 
inflationary energy scale to be A = 1.3 x 10 14 GeV in this 
model. One finds that these values imply 



-4.2 



A = lO 02 , g = 10 



-3.7 



(82) 



By recomputing the excluded region in the X-g plane for 
this value of v/M p , we find this point to be solidly ex- 
cluded. 

To get out of the excluded region, one should try to 
achieve smaller values of g or larger values of v/M p . One 
can show from the above relations that 



M„ 



ir 6 g s 



1/4 



m7, 



(83) 



Since A is fixed by the COBE normalization, lowering 
the string coupling is the only way to increase v/M p . 
Decreasing g s also has the desired effect of reducing g. 
The other parameter we can adjust is the warp factor, 
but there is not much leeway for decreasing a., because 
the density contrast in this model is given by [5Sj | 



1/3 



4/3 



(84) 



where C\ is a constant, and N e = 60. (Lower values of 
N e would require lower values of A, which push v down, 
contrary to what we want.) T3 cannot be pushed above 
the Planck scale, so a, can only be decreased by at most a 
factor of 10~ 3 / 4 . Although it also goes against the valid- 
ity of the low-energy effective supergravity theory upon 
which this model is based, in the interest of being conser- 
vative when deriving constraints, we will consider values 
of the warp factor which are this small. By scanning over 
g s and <Zj, we find the allowed and excluded regions in the 
plane of the warped string scale and the string coupling, 
shown in figure (|14|) . One consequence is an upper bound 
on the string coupling, 



9s < 10" 



(85) 



which is stringent, and smaller than normally expected. 
Our bound differs from those found in [l5j |. which also 
considered the nongaussian perturbations produced by 
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FIG. 14: Allowed and excluded values of the warped string 
scale and string coupling in the KKLMMT model, due to 
tachyonic preheating. The fiducial point corresponding to 
18211 is at (—1,-3.9), off the scale of this graph, but in the 
excluded region. 



tachyonic preheating. Their allowed regions are also ex- 
pressed in the plane of M s and g s , but differ markedly 
from ours. 

A caveat to this result is the potential for complica- 
tions due to the peculiar form of the tachyon action (|76() 8 . 
The dynamics of the open string tachyon action is differ- 
ent from the usual Klein-Gordon action which our anal- 
ysis assumes. In the case of <|76|) the potential is mini- 
mized at T = 00 and it takes an infinite amount of time 
for the condensation to complete (and the quantum me- 
chanical decay of the brane-antibrane system into closed 
string states invalidates this classical description on a 
time scale comparable to (ajMs)" 1 ). Thus the criterion 
((5a) 2 ) ~ v 2 defines the limit of applicability of our anal- 
ysis, rather than the end of the symmetry breaking pro- 
cess. Moreover, the field theory (|76[) leads to finite-time 
divergences in the spatial derivatives d{T near the core 
of the topological defects where the field stays pinned at 
T = [62, |65| , as well as finite-time divergences in the 
second derivatives didfT in regions where V — > 0, called 
caustics [6(| . It is not clear how these peculiarities might 
modify our estimate for Jnl, but we plan to investigate 
this in future work. 



8 The infiaton is also described by a DBI action (though the po- 
tential does not appear as a prefactor), however, slow roll ensures 
that it will always be a good approximation to expand the square 
root to obtain a standard kinetic term. In the case of the tachyon 
this may not always be a good approximation, particularly dur- 
ing the instability phase. 
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VIII. CONCLUSIONS 

We have carefully calculated the contribution from 
tachyonic preheating to the spectrum and bispectrum of 
the curvature perturbation in a model of hybrid inflation, 
and derived constraints on the model parameters using 
experimental limits on the power spectrum and nongaus- 
sianity in the CMB. The tachyonic contribution to the 
spectrum, which we parametrized by /l, has spectral in- 
dex n = 4, and the experimental limit on fa turns out to 
always be a stronger constraint than that coming from 
nongaussianity. Thus we do not expect to see nongaus- 
sian signals in future data sets, if hybrid inflation is the 
underlying model. Rather, one should look for an ex- 
cess of power at small scales to see the effects we have 
investigated. 

We have studied a simplified model of hybrid inflation, 
in which the tachyon is a real scalar field. This model 
is not viable as it stands, because domain walls will be 
copiously produced at the end of inflation, and they will 
quickly dominate the universe. We need to assume that 
there are some small operators in the Lagrangian which 
violate the symmetry a — > —a, and make the domain 
walls unstable. Although our results are strictly valid 
only for hybrid inflation with real fields, we have made a 
preliminary attempt to apply them to the warped brane- 
antibranc inflation model, by matching the coefficients 
of the lowest dimension operators in the corresponding 
Lagrangians. To the extent that this is valid, we find an 
interesting limit on the string coupling, g s < 1CP 45 . 

One achievement of this work is our demonstration 
that nonlocal contributions to the second order curvature 
perturbation, as defined herein, vanish when care is taken 
to consistently include all terms of a given order in the 
slow-roll expansion. The nonlocal terms which appear at 
intermediate steps when using the generalized longitudi- 
nal gauge have been the bane of some previous attempts 
to do similar calculations. They should not appear in 
any physical quantities since they violate causality. Any 
estimates made using the nonlocal terms are likely to 
grossly overestimate perturbations on large scales, since 
these terms are singular at small wave numbers. 

Some readers may feel uneasy about an event at the 
end of inflation being able to affect fluctuations whose 
wavelength corresponds to the horizon at the b eginning 
of inflation. However, as has been emphasized in |2a . l30| . 
there is no violation of causality here. Consider a collec- 
tion of comoving observers who are within the horizon at 
the beginning of inflation, but are about to lose causal 
contact with each other. If they agree to set up pertur- 
bations with small wavelengths in their Hubble patches 
at the end of inflation, this will induce large-scale cor- 
relations between different Hubble patches. Since the 
tachyon field is synchronized across these different Hub- 
ble patches it is able to amplify the large-scale perturba- 
tions in a coherent way. 

One avenue for future study is the more realistic class 
of models where a is complex and the defects which form 



are cosmic strings instead of domain walls. However, this 
is a much more technically difficult model to study than 
the case of the real scalar. In the latter case, the VEV 
of (7 remains zero even during the tachyonic preheating 
phase, because the field rolls in opposite directions on 
either side of each domain wall, and the VEV averages 
to zero. In the complex case, the modulus of er grows 
uniformly during preheating, and the randomly varying 
phase is only a gauge degree of freedom. We can no longer 
decouple the equations for the fluctuations of the metric 
and the inflationary fields when a has a VEV. It thus 
remains an interesting question whether the realistic hy- 
brid inflation model also gets significant nongaussianity 
from tachyonic preheating |67|. 
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APPENDIX A: Generic Field in deSitter Space 

In this appendix we review the solutions of the Klein- 
Gordon equation in deSitter space 0, |3|J . If we consider 
the limit of deSitter space e = rj = then ipo,H are 
constant and both the tachyon and the inflaton first order 
fluctuations satisfy an equation of the form 



6 + 3H£ k + 



r,2 



6 = o 



(A-l) 



or, in terms of conformal time 

^'l + 2H& + [k 2 + a 2 m\] & = 



where a = e Ht 
fk = a£,k the equation simplifies to 

n 2 



-1/Ht and H = — 1/r. Defining 



ft 



1 



(A-2) 



This is the variable in terms of which the action is canon- 
ically normalized. Choosing the Bunch-Davies vacuum 
corresponds to normalizing the solutions so that [39| 



fk(r) 



-ikr 



/V2k 



(A-3) 



on small scales k ^S> aH or equivalently — kr ^> 1. The 
general solutions to (|A-2|) are 



fk(r) 



Cl 



(fc)i?W(-fcT)+ C2 (fc)ij( 2 )(-fcr) 



(A-4) 

w here the ord er of the Hankel functions is v — 
9/4— m^/H 2 . This analysis makes no assumptions 
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about the size of m^/H and all the formulae we derive 
are valid for arbitrary complex v unless otherwise stated. 
Notice that — r > for all t so that the arguments of 
the Hankel functions are always positive definite. We re- 
cover the desired asymptotics (|A-3ll on — fcr <C 1 with 
the choice 



ci(fc) 



exp 



The solution becomes 



fk(T) = — - — exp 



in f 1 
— \v+- 
2 V 2 



ca(ft) = 0. 



H^i-kr) (A-5) 



or, in terms of cosmic time and the original field variable 



&(*) = o 



2 V a 3 H 



exp 



?.7T 
Y 



1 

" + 2 



' aH 



(A-6) 

To get an intuitive feel for this initial condition we con- 
sider the asymptotic behaviour of l|A-5|l on large scales: 



fk(r) 



y/ — TTT 



■ exp 



ITT 

Y 



i 

V+ 2 



(A-7) 



1 



kr 

Y 



kr 

Y 



as — kr — > 0. To simplify further we should focus on 
either m^/H > 3/2 or m^/H < 3/2. 

We consider a field which is heavy during inflation 
7Ti£ 3> H so that i/ w im^/H is pure imaginary. How 
do the two terms in the square braces of ljA-71) com- 
pare? The functions {—kT) imi ' H are oscillatory so the 
relative magnitude of these two terms depends only on 
the gamma function prefactors. Using the results (for (3 
real and arbitrary complex z) 



T(l + z) = zT(z), |r(l+i/3)| = 



7T/3 



sinh(7r/3) 



from the theory of gamma functions we find that for 
me/H > 1 



r(l + im$/H) 
\r(im £/ H)\ 



'-^-expf^i) 
27rm 5 V 2H J 



l2irH 



■ exp 



so that the first term in the square braces on the second 
line of l|A-7|l dominates at large ra^/H. Going back to 
the original variable we find 



on large scales. 



-3/2 _ 



2 3 / 2 , 



for rri£ 3> i? 



We consider now a field which is light during inflation 
77i£ <C H. In this case f « 3/2 is pure real and the 
first term in the square braces in l|A-7(l goes to zero. We 
recover the scale invariant spectrum 

|&| ~ -^L= for m 4 < H (A-8) 



on large scales, as expected. 

Notice that for any value of /H we have the small 
scale asymptotics 



2k 

by constuction. 

During inflation a(t) grows exponentially so that the 
large scale fluctuations of any massive field are damped 
as a~ 3 / 2 while the large scale fluctuations of any light 
field are approximately constant. On small scales the 
fluctuations of all fields get damped as a -1 . 



APPENDIX B: Perturbed Einstein Equations and 
the Master Equation 

Using Maple we have carefully verified the results 
of |32| for the perturbed Einstein equations and the 
master equation. Here we briefly review those results 
relevent for the computation of (( 2 \ We present only 
the S^G° Q = k 2 5^T°, diS^Gi = kV^T* and 
SjS^Gl = K^SjS^T? equations since the second order 
vector and tensor fluctuations decouple from this system. 
In the case that <jq = the second order tachyon flucta- 
tion 5^ a decouples from the inflaton and gravitational 
fluctuations. Analogously to the first order fluctuations, 
the Klein-Gordon equation for the inflaton fluctation at 
second order S^tp is not necessary to close the system 
of equations. In the subsequent text we sometimes in- 
sert the slow roll parameters e and r\ explicitly though 
we make no assumption that they are small. We also 
introduce the shorthand notation m 2 , = d 2 V/d<p 2 and 
m 2 = d 2 V/da 2 and assume that d 2 V/dadip = 0. 



The second order (0, 0) equation is 
mip' (2) + (3 - e)H 2 <t> {2) - 



~2 r 



Ti 



(B-l) 



where Ti is constructed entirely from first order quanti- 
ties. Dividing Ti into inflaton and tachyon contributions 
we have 

Ti = Tf + TI 

where 

Tf = 4(3 -e)H 2 (V 1} ) 2 + 2 K Vo^ (1) <5 ( V 

+ ^d k d k ^ +3^y + 3(v<t>^y (B-2) 
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and 



Tf = 



Y 



(B-3) 



Notice that we have decoupled ip^ and from the 
inflaton perturbation S^tp. It remains now to express 
-0( 2 ' in terms of 0^ 2 '. To this end we subtract the in- 
verse laplacian of (|B-4|I from (|B-7|I and again use (|B-4|) 
to eliminate 5^(p which gives 



d k d k 



The divergence of the second order (0, i) equation is 

d k d k [<// (2) + n^] = ^d k d k 8^ V + T 2 (B-4) 

where T2 = T 2 + T 2 is constructed entirely from first 
order quantities. The inflaton part is 

T£ = 2n 2 p , d l ^d l 8^^j + k 2 ^ (,5 ( V 3M (1 V) 

- 89, (^W^) - 2ft (V 1 )^ 1 ') (B-5) 

and the tachyon part is 

T a 2 = n 2 d l (6 il) a'd i 5( 1) a S ) . (B-6) 

The trace of the second order equation is 



0( 2 ) - </> (2) 



= T 3 - 3A _1 T; - mA~ l T 2 



or, equivalently 



(2) 



A _1 7. 



Following the notation of [32] we have defined 



(B-ll) 



(B-12) 



/(2) 



which can be split into inflaton and tachyon components 
7 = 7 V + 7 CT in an obvious fashion. Our T 2 is related to 
the quantities a, (3 defined in [32j by 

T2 = n 2 (i — a. 

Now, using (jB-llfl to eliminate from ljB-10|l gives 
the master equation 

(p"^ + 2( V - e)H^ 2) + [2(77 - 2e)H 2 - d k d k ] ^ 
= T 1 + A- 1 T 2 -2(2 + e-? ? )WA- 1 T 2 -7 
- (l + 2e- 277)7iA- 1 7 ' + A~V 



3H^' (2) +3(3-e)H 2 (2) 

3^ 

2 



(B-7) 



where T3 = T3 + T3 is constructed entirely from first 
order quantities. The inflaton part is 

2 



Inserting explicitly the expression for 7 (|B-12() this can 
be written as 

<P"^ + 2(7? - e)UP' {2) + [2(77 - 2e)H 2 - 8 k d k ] ^ 2) = J 
where the source is 

J = T 1 -T 3 + 4A- 1 T' 2 + 2(l-e + T ] )HA- 1 T 



1 3 



12| 

3k 2 



(3- e )W 2 (^)) - 6«Vo0 W * ( V 

a 2 m 2 (^V)'-^(^V 
+ 3 (>«) 2 + S4> (1) d k d k 4> (1) + 24W0<V (1) 

+ r(v^) 2 (b- 

and the tachyon part is 
3 



-3 "I--"- 1 J- 2 

A- 1 7"-(l + 2e-27 7 )WA- 1 7'. 



2 

(B-13) 



2 1 



2„.2 

-a m„ 



(B-9) 



We now proceed to derive the master equation. Adding 
(|B-1|) to the inverse laplacian of the time derivative of 
(IB-4|) and then using (|B-4|) to eliminate S^tp yields 

y/' (2) - (1 + 2e - 2r 1 )H^ 2) + H<t> l{2) 

- 2(2e - V )H 2 (j} (2) - d k d k ^ (2) = Ti + A" 1 ^ 

- 2(2 + e-r/)WA- 1 T 2 . (B-10) 



We can split the source into tachyon and inflaton contri- 
2 butions J = J v + J a in the obvious manner, by taking 
the tachyon and inflaton parts of Ti, T 2 , T3,7. 

We now derive some results concering the tachyon 
source terms which will be useful in the text. First we 
consider the tachyon contribution to 7 l|B-12|l : 

^ 7 CT = Y3 — 3A _1 9 T T 2 — 67YA _1 T 2 . 

Using equations 1|B-6(I and i|B-9() we can write this as 

- iK 2 A- 1 d T d l {5 {1) a'd l 5 {1) a) 

- G^HA^d^a'd'S^a) (B-14) 
We now write 7^ as 



^d k d k (5^a>) 2 -\d k d k (d t 6^od*5^o) 



^a 2 m 2 a d k d k (S^a) 2 - 38^(8™ a' 8*8^ a) 
6Hd t {8 (1) a'd l 5 (1) a) 
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and, after some algebra, we have 



la 



7<7 



The last two lines can be simplified using the equation of 
motion for the tachyon fluctuation 

# V + 2HS^a + a 2 m 2 J^a = d k d k S^a 

which gives 

- 3d l (d k d k S (1) a)d l 8^a - ?,d l d l 5 {1) ad k d k 8 w a 

\d k d k (ft^W^V) 

- 3d l i K d k d k S w ad l S {1) a^ (B-15) 

This result has also been derived in ■ We now consider 
the tachyon contribution to the source. We take 

j* = Tnj r -T£ + 4A- 1 d r T£ + 2(l-e + ?7)?iA- 1 T£ 
+ A- 1 7 "-(l + 2e-27 ? )HA- 1 7 ; 



and, using (|B-3|) . IB-6|) and (IB-9(I . we have 

+ 2 K 2 H{1 + f]- e)A~ 1 d i (6^a'd l 5 {1) a^ 
- H{l + 2e-2r 1 )A- 1 1 ' a + A- 1 ^. (B-16) 



APPENDIX C: The Inflaton Contribution to < (2) 

In this appendix we consider the calculation of the in- 
flaton part of the second order curvature perturbation 

(2) 

Q> using results from appendix B and following closely 
the calculation of the tachyon part of the second order 
curvature perturbation which is discussed in subsection 
HVDl 

1 1 (2) 

The first step to computing Q, is to use the Green 

function 14411 to eliminate the dependence of (}p on </>( 2 ' 
(from the first line of l|48l) which is equivalent to \b'6\i ) in 
favour of (j)^, S^ip. This procedure follows closely the 
analysis of subsection IIV Dl and in particular (|56|) still 
holds though here we are interested only in the inflaton 
part of the soucre J v . 

Having removed the explicit dependence of £( 2 ) on the 
second order metric perturbation tp^ we next eliminate 



</>W, S^tp in favour of Using the solution i|22|l and 
the first order constraint equation l|15|) one may verify 
that on large scales 



c 



(1) 



n e 



and 



.^0,(1) 



On large scales the first order curvature perturbation is 
approximately constant since |42( 



'(i) _ 



using the fact that oo = (so that there are no 
anisotropic stresses, whose absence guarantees the con- 
servation of C*- 1 ) on super-Hubble scales). Thus on large 
scales we have 

= -(26-rMC (1) - 

It is straightforward to compute the last three lines of 
gEJ using the fact that Q { J } « -<^C (1) /H and Q'j 1] « 
— (2e — ??)¥>oC on large scales. The result is 



C (2)g 





fam v \ 


l 


V H ) 



2 + 2e - 2-q 



(C W)2 



(C-l) 



where we have dropped terms which are higher order in 
slow roll parameters or which contain gradients. Notice 
that the quantity (a m v /TC) 2 is first order in the slow roll 
expansion because 

1 1 d 2 V 

k 2 V df 2 



1 



1 2 

1 m-tp 
k 2 V 



1 / a m v s 2 



To (IC-lll we must add the contribution coming from the 
fist line of (|48(1 which can be written explicitly in terms of 
first order quantities using equations H53fl and H56[) . Here 
we consider only the particular solution for (f>^ due to 
the inflaton source J v . 

In order to compute J v we first study the quantities 
Ti, T2, T3, T4, 7 which are defined in appendix B. On 
large scales and to leading order in slow roll we have 



yip 
L 1 



12e 2 



(- 



n 



n 2 (c {1) ) 2 



The quantity A can be written as (see equation 39 

of |a2) 



A -i T 



Xyip 



^(e~ v )H(5^^) 2 



m(<b {1) ) 2 



(i)^'(i) 



—A' 1 (d i d k d k (t> (1) d i S^ip 



+ d k d k (j> (1) d l d l 5 {1) <p) 
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The inverse laplacians on the last two lines contribute 
only to the momentum dependence of f% L which we ne- 
glect. On large scales and in the slow roll limit we have 



A" 



1 2 



(4 e 2 -^)H(C (1) ) 2 



where the • • • denotes momentum dependent terms. On 
large scales and in the slow roll limit we also have 



1 3 



36e 2 - 3e 



/am v \ 
V U ) 



h2 



The quantity T 4 defined in l|54|) depends only on the 
inflaton fluctuation and is given by 



T 4 = T£ « (6e 2 



2^)^(C (1) ) 2 



on large scales and in the slow roll limit. Using these 
results one may readily verify that j^/Tl 2 is third order 
in slow roll parameters 

7 ^0( £ 3 )H 2 (C«) 2 . 

Using these results one may verify that the only term 
in (|53|1 which contributes at lowest order in slow roll pa- 
rameters is the term proportional to /eli 2 . Thus, the 
contribution to the second order curvature perturbation 
due to the first line of (|48|l is 



c 



(2) 



4e- 



1 / a m v y 
3 V H J 



(C (D 



3 v n 

Adding all the contributions together we find 

The contribution 2(Q 1 ') 2 stems from using the Malik and 
Wands ^lj definition of the second order curvature per- 
turbation. It can be related to the definition of Lyth and 
Rodriguez (which also agrees with Maldacena 0) 
using 



2(C (1) ) 2 . 



The Lyth-Rodriguez curvature perturbation, due to the 
inflaton up to second order, can thus be written as 



LR 



c 



(1) 



/M (1) ) 2 



where 



/£l=«(2^-6 £ ). 



In writing C t f jR above we have suppressed the homoge- 
neous k = mode of ( which should be subtracted to 
ensure tha t (( ) = 0. This result differs from previous 
studies E5 by a factor of two. The calculation of 
E3 takes into account the effect of nonlinear evolu- 
tion aswell as the effect of computing the bispectrum in 



the vacuum of the interacting theory, as opposed to the 
vacuum of the free theory. Our calculation does not con- 
sider the effect of the change in vacuum which is the same 
order of magnitude. Thus we should not expect to repro- 
duce exactly the results of 0] • The change in vacuum 
will not change the calculation of the tachyon part of the 
curvature perturbation since £W does not depend on 
S^'a and hence contributions to £ CT due to the change in 
vacuum will be higher than second order in perturbation 
theory. 



APPENDIX D: Fourier Transforms, Mode 
Functions and Inverse Laplacians 

We define the Fourier transform of some first order 
quantity 5f(t,x) by 

d 3 k 



6 f(t, 2) 



(27r) 3 / 2 
d 3 x 

(27t) 3 / 2 ' 



c 6f k (t) 



(D-l) 



ik-.r 



x 5f{t,x) 



where the Hermiticity of Sf(t,x) implies that 5fk(t) 
8f-k{ty so that we can define 



5fk(t) = «*&(*) + aL$_*(t)* 



(D-2) 



where a k is an operator and Sfk(t) is a c-number valued 
mode function. We then re-write the Fourier transform 



as 



5f(t,x) 



d 3 k 
(2tt) 3 / 2 

d 3 k 
(2tt) 3 / 2 



In this form it is clear that a k is the usual creation oper- 
ator satisfying 



and one should expect 
Mr) 



5 {3 \k-k') 



-ikr 



2k 



on small scales, which corresponds to the Bunch-Davies 
vacuum a k \0} = 0. This is consistent with the usual 
definition of the power spectrum in terms of the two- 
point function 



(0\6f(t,x)\0) 



so that 



Vj{k) 



2^2 



d 3 k 
dk , s 



l&(*)| S 



24 



We now discuss the Fourier transform of the tachyon 
source. Typical terms in the source have the form 

J«(i,f) = b(t)5f(t,x)6g(t,Z), 

jW(t,x) = 6(t) A" 1 a)<fo(*.2)]> 
jW(t,x) - b(t)A- 2 [6f(t,x)6g(t,x)] 

where 5f, 5g are some first order quantities and b is con- 
structed from zeroth order quantities. The Green's func- 
tion for the laplacian is defined as appropriate in the 
absence of boundary surfaces 

(A-\f)(t,x)= J d 3 x'G(x-x')f(t,x") 

where 

G(X-X') : 

In Fourier space we have 



1 



1 



4tt \x — x! 



G{x,x') = 



G k (x') 



d 3 k 

(27T) 3 / 2 ' 



c G k (x') 



-ikx' 



(2tt) 3 / 2 k 2 ' 
In Fourier space the souce terms contain convolutions 

4 1] (t) = b(t)(6f*5~g) k (t), 



4 2) (t) = -b(t)^(6f*5g) k (t), 

4 3) (*) = Kt)±(5f*6g) k (t). 

where f k , g k are operator valued Fourier transforms de- 
fined as in l|D-l|) . These may be related to the mode 
functions as (|D-2|I . Finally, we have defined convolution 
by 



1 



(Sf*8g)M) = 



d 3 k' 
(27r) 3 / 2 



skm k -p(t) 



APPENDIX E: Construction of The Tachyon 
Curvature Perturbation 



where the second order Sasaki-Mukhanov variable is 



Q< 2) - <5 (2 V + f^ (2) 

+ inflaton contributions 



(E-2) 



Inserting (|E-2|) into IjE-lfl and using the (0, 0) Einstein 
equation (|_B-lfl to eliminate the contribution ip^S^tp' + 
a 2 dV/d(pS^ 2 ' ip gives 



c?> = - 



cH 



i 



d k d k 



7a 



A- 



7<7 



3-e 
1 1 
3-eeH 



eH 2 
2 1. (E-3) 



where we have also eliminated ip^ in favour of cj)^ 
using (|B-11|) . Notice that using (|B-1I) introduces a 
term proportional to to the curvature perturbation 
which cancels the contribution proportional to (5^a') 2 + 
a 2 m, 2 (S^ a) 2 on the second line of (|E-1|) up to a gra- 
dient term (VS^a) 2 which can be neglected on large 
scales. The result <|E-3I) is valid only on large scales but 
we have not yet assumed slow roll. In (|E-3|I it is un- 
derstood that denotes only the particular solution 
due to the tachyon source, J" . We now use the Green 
function l)44|) to solve for in terms of the tachyon 
source (IB-16|) . We work only to leading order in slow roll 
parameters. We also work in the large scale limit. To 
lowest order in e, r\ and up to order k 2 we can write 



eH 



(2) 



k^ k 
3eH 2 



dr'e(r-r')J k (r r ) 



6 



5 o » 2 i 
— t t — r 

6 3 



(E-4) 



using the Green function 1441) and the relation H45fl . 
Equation i|E-4(l allows us to eliminate the dependence 
of (i 2 ^ on <j)^> . At leading order in slow roll parameters 
the tachyon source is (see (|Ij-13|) ') 



In this appendix we include technical details about the 

(2) 

contraction of Q- using the Green function for the mas- 
ter equation. For clarity we repeat some details which 
are included in the text. We consider only contributions 
to the equations which depend on the tachyon field and 
we remind the reader that is independent of a. The 
second order curvature perturbation is (see (148151)1) ) 



c 



(2) _ 



1 



1 



3 - e (^) 2 

1 1 



M 2) 



3 - e (^) 2 
inflaton contributions 



S^a'Y +a 2 m 2 



(>vy 



(E-l) 



J" 



+ a-^ 7ct 



4A _1 a T T5 



2WA~ 1 Y? 



(E-5) 



We must now insert (|E-5|) into (|E-4|) . perform numer- 
ous integrations by parts, and then insert this result into 
(|E-3|) . We evaluate term-by-term the last line of l|E-4|) . 
The first contribution (the term proportional to A; ) is 



*'i 2) 



2J,(2) 



eH e 3eH 2 ' e 
Inserting (|E-5|) . noting that H(t') 



dr 



H(r') ' 
T/t' at leading 
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order in slow roll and integrating by parts gives 



eH 



dr' 



ZeU 2 



1 

9 

e 



W(r') H(t' 



W(t') 

6A- X T£ 



H 



H 



(E-6) 



where the terms under the dr' integral are evaluated at 
t' while the terms in the square braces on the second 
line are evaluated at r. The • • • denotes constant terms 
evaluated at r = Tj which arise from the integration by 
parts. Since our interest is in the preheating phase during 
which the fluctuations S^'cr are amplified exponentially 
we can safely drop these constant terms. 

The second contribution to the last line of (|E-4(I has 
the form 



A 2 ) 



3eH 2 



1 

3 6e 



k 2 J° 
dr' * Jk 



H(t>) 



f\3 ■ 



In evaluating this we need only consider terms in k JZ 



which are not suppressed on large scales. Using (|E-5|) 
one may verify that 



k 2 r k = -i'' h 



IT, A' 



(E-7) 



on large scales (recall that Yf is a gradient, see p}-60 ) so 
that we have 



A 2 ) 




K 


_t 


eH 










dr' 






i 





3eH 2 




2T? 

3 +6A - 

3 H 


1 2 


1 ^ +2 A 


— 1 "y~cr " 
1 2 


3eH 2 + 


n 



dr 



k 2 T 



2 7<r 

k 



H(t>) 



l\3 



(E- 



Equation l|E-7|l shows why it was necessary to include 
the k 2 terms in the large scale expansion of the Green 
function l|E-4() . noting that 7 CT may be written as (|B-12() 
and comparing to (|E-5|I we see that k 2 Jk contains terms 
which are of the same size as those in J/~, on large scales. 

(2) 

Thus a consistent large scale expansion of Q requires 
that we work up to order k 2 in the expansion of the Green 
function. 

The third contribution to the last line of (IE- 41) is 



eH 

5 

OV 



3eH 2 



5 
6e 



dr' 



k 2 JZ 



H 3 



H{t')H{t) 2 
(E-9) 



using the same procedure as above. 



Finally, the fourth contribution to the last line of ljE-41) 



is 



A 2 ) 

eH 
2 



4 2) 



fc 2 4 2) 

3eH 2 



lie 



dr 



H{rf 



3eH( 



A^TS 

) £ 

' H 



3 H 2 3 H 3 



(E-10) 



Summing up l|E-6jl, (|E-8|) . (|E-9|) and (|E-10|) and inserting 
the result into (|E-3|) gives 



(2) ~ 



dr' 



1 1 

~H{t') T 3W(t') 



1 Tf 



c 



Now, using equations l|B-3(l and l|B-9j) we can write this 
in terms of the tachyon fluctuation as 



znjr') 2 • 

3 H{t) 3 3 



& - ^ f dr' 

e J-l/a t H 

- a 2 ml(5^a) 2 



{5^a'f H{r'f 



H{r>) H(t) 



(*V) a 

(E-ll) 



APPENDIX F: The Adiabatic Approximation 

Here we consider the construction of curvature per- 
turbation in the case in which the tachyon mass varies 
adiabatically and demonstrate explicitly that the result- 
ing bispectrum is scale invariant. If the inflaton rolls 
slowly enough then the tachyon mass will evolve in time 
extremely slowly and will remain close to zero for a sig- 
nificant number of e-foldings. We define the slow-roll 
parameter for the a field as 



Xv 



M„ 



(F-l) 



For example, if ~ 30, then [rj^l has the same value at 
the beginning and end of inflation. We will denote the 
magnitude of the value of r\ a at the end of inflation (and 
preheating), t = i/, by 



Vf = \Va(tf)\ 



(F-2) 



The condition to have a slowly rolling tachyon during 
2N r e-foldings of inflation is therefore 



T)f < I 



< 



Xv 6 



32-/V*M* 



(F-3) 



Combining this with the condition (|24|l . it follows that 
the coupling g must be extremely small if >, I, 



< 



io- 



v 



(F-4) 



v 
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Note that a small coupling g does not require fine tuning 
in the technical sense, since g 2 is only multiplicatively 
renormalized: (3{g 2 ) ~ 0(g 2 X, 5 4 )/(16tt 2 ). That is, if g 
is small at tree level, loop corrections will not change its 
effective value significantly. The bound l|F-4(l insures that 
the tachyon will remain light compared to the Hubble 
scale during N* e-foldings of inflation. 

The combination of (|F-4|) and 1|27[) . eliminating g, gives 
an upper limit the tachyonic mass scale required for the 
consistency of our assumptions, 



Xv 2 < 



4 x 1(T 

ivT 



■ Mi 



(F-5) 



We now turn to the adiabatic approximation for the 
tachyon mode functions. In this limit the tachyon mass 
is varying adiabatically, so we can determine the behavior 
of the fluctuations from the deSitter space solutions by 
replacing the dependence on a constant tachyon mass 
with time- varying one. The deSitter solutions for massive 
fields are given in appendix A. We are interested in the 
long-wavelength growing modes, which behave like 



8 {1) a{x) = 



d 3 k 



H 



(2tt) 3 / 2 



■X-hry 



'e lkx a k 



h.c. (F-6) 



in conformal time. Since r\ a is slowly varying, we can 
replace it with its time-dependent value (|F-1|) . that is 



M 2 

t? ct (t) = 877 — In 

1 T-^ 



a(r) 



(F-7) 



Note that ln(et*/a(Tf)) = — N* at the end of inflation, 
so rjf — 8r]N*M 2 /v 2 . We can use lj5U|) . which is the cri- 
terion for the end of tachyonic preheating and inflation, 
to determine iV» in terms of the parameters of the hy- 
brid inflation model. Using l|F-6() . the dispersion can be 
computed at the end of inflation, r = Tf 



(S w a 2 (x)} 
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where rjf — \T)o-(T~f)\ > 0. The infrared and ultraviolet 
cutoffs are defined in the same way as is needed to make 
the dispersion of the scale-invariant inflaton fluctuations 
finite: 9 \koTi\ = 1 at horizon crossing of the largest scale 
modes visible today, so that |fcor/| = e~ N " with N e ~ 60 
being the total number of e-foldings of the visible part of 
inflation. Similarly |fc/7y| = 1 at the end of inflation: kf 
represents the last mode to cross the horizon before the 
end of inflation. Even though the integral converges in 



9 The first order inflaton contribution to the curvature perturba- 
tion's two-point function has the form ((Q, ,(x)) 2 ) ~ J d 3 fe/fc 3 , 
and also requires cutoffs 



the absence of the cutoff, later on we will need to know 
the correct behavior of (|F-8I) when 77/ is small, which is 
not correctly represented by the limit fcmax — * 00. 

Eq. HF-8f) can only describe the growth of the tachy- 
onic fluctuations until eq. (|30|l is fulfilled, and the per- 
turbative description breaks down. This signals the end 
of inflation, and allows us to determine the number of 
c- foldings during preheating, N», in terms of the funda- 
mental parameters of the theory. Using (|F-8|) and ijF-lfl . 
the implicit relation 



Xv 6 



64M£m 2 



In 



1 + 768tt 2 AT* 



M p b m 2 



(F-9) 



follows. Using the COBE normalization (|24(1 to eliminate 



niip , this can be rewritten in the form 
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This can be solved iteratively on the computer; using an 
arbitrary guess for and repeatedly evaluating l|F-10(l 
gives a rapidly converging answer. 

For consistency, we should check whether we have the 
freedom to make TV* > 1 given other restrictions on the 
parameters. From (|27J), gVX > 230(v r Xv/M p ) 3 , and tak- 
ing N e ~ 60, we get an upper limit on iV*, 



7 M l , 

1 X~^ 



0.05 IniV* 



(F-ll) 



This would pose a problem for our scenario if it was 
not possible to take Xv 2 to be sufficiently small, say 
Xv 2 /M 2 < 10 -8 . However such a value is consistent with 
the restriction l|F-5J) . so there is no new constraint com- 
ing from (|F-8J| . This demonstrates that it is consistent 
to assume that the tachyon remains light until the end of 
inflation, provided that g is sufficiently small, eq. (|F-4|) . 

We now perform the time integral in the result 157() 
using the knowledge of how S^a behaves in the long- 
wavelength limit which will be of interest for cosmologi- 
cal observables, l|F-6|l . According to this result, S^a' ~ 



I tJ^ tj] whereas 



^mi/(H 2 T 2 ) =3r? CT . There- 



fore, since we are working in the regime where the 
tachyon is light, r\ a < 1, the term with m 2 in l|57|l dom- 
inates over the terms with time derivatives. Let Tf be 
the final time in the upper limit of integration in 1|57|) . 
and r* the time at which the tachyonic instability be- 
gins. For concreteness, we will at first assume that there 
are equal numbers of e-foldings (N* ~ 30) before and 
after this time, so that the slow-roll parameter for the 
tachyon is time-antisymmetric, ^(Tj) = — ^(t/), since 
this is the way of having a light tachyon during inflation 
which minimizes the tuning of the tachyon mass. Chang- 
ing variables to x — — ln(r'/r*), and writing r\ a = —dx 
with d = rjM 2 /v 2 , eq. I|57|l in Fourier space takes the 
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form 



e\k) ~ 3d^ J (S^a) p (8^a) k . p (F-12) 



dx \ P n\- dx \{k ~p) T ,\- dx e 2dx2+3 ^- N ^ + 0(d 2 ) 

N, 

where S^a p is the part of 5^o p left after dividing 
by the only time-dependent factor, (— pr)''". In fact 
H/ (2p 3 ) 1 ^ 2 is the perfectly scale- invariant part of 
the perturbation, as far as its contribution to the power 
spectra of fluctuations, while the (-pr) 11 " gives rise to 
the small departures from scale invariance. 

The integral (|F-12[1 is dominated by contributions near 
the final time, so one expects it to behave roughly like 
the value of the integrand at x = iV*. By numerically 
evaluating (|F-12|) and studying its dependences on , d 
and fcr* , we are able to deduce the following semi-analytic 
fit: 



C 



(2) 



1_ r 2y f N, 



1 



\pr* 



(2tt) 3 / 2 

Vf 



(6^a) p (S^a) k _ p 
(F-13) 



where r\f — \f)a(jf)\ is the magnitude of the final value of 
r\ a . Notice that \pr„\ = p/(a^,H), so that a mode which 
crossed the horizon at the beginning of inflation (satis- 
fying \pTi\ — 1) would have \pT*\ = e~ N ". The second 
order perturbation at these scales therefore grows by a 
factor of e ir)fN ' relative to scale- invariant perturbations, 
assuming \k — p\ ~ p. This is the growth which is due to 
the tachyonic instability. We remind the reader that this 
growth saturates and our perturbative treatment breaks 
down when S^a ~ v/2, whereas initially S^'cr ~ H, so 
the result (|F-13(I is only valid when e 471fN " <, v/H. 

We can also consider a different case, when horizon 
crossing occurs at the time r* , and then N* ~ 60 instead 
of 30. We find that the formula (|F-13|1 is essentially un- 
modified, because of the fact that the integral is domi- 
nated by its late time behavior, so that the contributions 
from the period r < r* are negligible. In this case, the 
total growth in Q 2 ' comes from the factor e 2l7/JV * and not 
from | fcr* | , since the latter is of order unity. Nevertheless 
the total amount of growth in the two cases is the same, 
~ (e 60l?/ ) 2 , and its range of applicability is limited by the 
same factor (v/H) 2 . 

We now calculate the leading contribution to the bis- 
pectrum of the second-order curvature perturbation (l(jl|) . 
Since there is no momentum dependence in the result 
for C (2) 



(|F-13(1 for Q ' , this three-point function is straightfor- 
ward to compute, using the massless free-field two-point 
functions 



H 2 



-.SWfa+Qi) (F-14) 



Carrying out the contractions of pairs of fields which 
contribute to the connected part of the bispectrum, one 



finds eight terms, which are compactly expressed in terms 
of the wave numbers pi and qi = ki — pf. 



&P1+P2 (^92+P3 ^93+91 + <^92+93^P3+9l / 



(F-15) 



+ <5pi +92(^2 +93 <^P3+9i + <5p2+P3^93+9i) + \Pi 1i) 



where 

f(jH,ki) =d*\n\- 2 ^\ Pl \- 3 / 2 -Vf \ kl ~ Pl \- 3 / 2 -Vf (F-16) 

and d t = H 2 K 2 rife 2r,fN */(2e). Each triplet of delta func- 
tions contains the factor 5k 1 +k 2 +k 3 expressing the trans- 
lational invariance of the bispectrum. The remaining two 
delta functions can be used to perform the integrals over 
Pi and p 2 , leaving the p%= p integral. 

The resulting bispectrum can be shown (with help from 
Maple) to reduce to the sum of two terms, 

B(ki) = MrA- evf 4 J ^\p\- 3 - 2 ^\p-h\- 3 ^ 

(F-17) 



(2tt) 

{\p + k 2 \- 3 - 2 ^ +\p + k 1 \- 3 - 2 ^) 



This expression is manifestly symmetric under k\ «-> k 2 , 
but it can also be shown to be symmetric under inter- 
change of any two k^s by an appropriate shift of the 
integration variable. The integral converges as p — > 00, 
but since r/f > by definition, it diverges with the small 
power l/p^f near p — 0, and we must introduce an in- 
frared cutoff, ko, representing a scale a few e-foldings 
beyond our present horizon. This procedure was also 
necessary for the dispersion of the fluctuations ((S^a) 2 ) 
in eq. (|F-8|I . But unlike in l|F-8|l . we must cut off the 
integral not just for \p\ < ko, but also for |p — ks\ < fco, 

+ < kg, and Ip + fel < k$. The physical reasoning is 
however the same: there should be no observable effects 
coming from fluctuations whose wavelength is far beyond 
our present horizon. 

By dimensional analysis, one can see that if all 
wavenumbers have the same magnitude, ki = kki, then 
the bispectrum goes like fc~ 6 ( 1 +''/)l' r *l f times a dimen- 
sionless function of the directions, ki ■ kj, of k/ko and 
of rjf. For simplicity, we evaluate it at the symmet- 
ric point fci = (-1/2,^3/2,0), k 2 = (-1/2,-^3/2,0), 
k-j = (0, 1, 0), where the three wave vectors form an equi- 
lateral triangle. The integral is evaluated numerically to 
obtain the result 



B(k t ) = p/a(fc/*d,»7/) 



H 2 K 2 r) f e 2r if N x " 



2ire\kTJ 2r if 



(F-18) 



We plot the function f 3 (k/k ,r]f) in figure fT31 which 
shows strong sensitivity to the presence of the cutoff for 
modes with k < ko, but weak scale-dependence for modes 
just a few e-foldings shorter in wavelength. Since observ- 
able quantities should not be sensitive to the exact value 
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of the cutoff, the region with ln(fc/fco) > 1 is the physi- FIG. 15: Log of fi,{k/ka, rjf) versus In (fc/fco) for several values 
cally relevant one. In this region, a reasonable fit to the of 77/, and the large-fc fit i|F-19(i . 
result for rjf <, 0.1 is given by 

(F-19) 

which is also plotted in figureEl It is clear, then, that in 
the case of an adiabatically varying tachyon mass there 
is only a small breaking of scale invariance. This break- 
ing of scale invariance can be neglected when evaluating 
the bispectrum near k = ko, corresponding to the largest 
presently observable scales. Our infrared cutoff is effec- 
tively at wavelengths a few e-foldings larger than this 
scale. 



6 




'0123 
In k / k„ 



We have computed the nonlinearity parameter /jvl 
using the approach detailed in this appendix and have 
found that large nongaussianity is produced for all values 
of the parameters which are consistent with the assump- 
tions we have made. 
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